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Abstract 


In  order  to  expand  shear  testing  to  different  nominal  strain  rates  a  new  experimental 
technique  iias  been  developed  in  the  Laboratory  of  Physics  and  Mechanics  of  Materials,  Metz, 
France  (lPMM).  A  modified  specimen  geometry  of  double  shear  has  been  designed  (MDS 
specimen  1.  'fhe  MDS  specimen  with  flat  faces  is  loaded  directly  by  a  bar  projectile  which  can  be 
accelerated  to  a  aesired  impact  velocity  by  a  gas  gun  (1  my's  <  Vo  ^  200  m/s),  where  Vq  is  the 
impact  velocity.  Suce  the  .specimen  is  backed  by  a  long  tube  which  transmits  the  force  imposed  on 
the  MDS  specimen  the  shear  stress  history  can  be  exactly  recorded  by  the  transmined  elastic  longi¬ 
tudinal  wave.  Tne  shear  displacement  of  the  MDS  specimen  is  measured  directly  by  an  optical 
sage.  The  experimental  setup  which  consists  of  the  air  gun,  transmitter  tube  (Hopkinson  mbe)  and 
ooucal  displacement  gage,  permits  for  a  wide  variation  of  the  nominal  strain  rates,  typically  10-  s'- 
<  r  <  10-  S'k  On  the  other  hand  there  is  no  limitation  of  maximum  shear  strain. 

Tne  MDS  soecimen  geomemy  has  been  calibrated  with  the  FEM  and  ".AB.ACUS"  software. 

Evaluation  of  dispersive  effects  of  elastic  longitudinal  waves  in  mbes  has  been  studied,  see 
.Appendix  n®  2  to  the  Repon.  A  cicsed-fonn  solution  has  been  obtained  for  a  semi-inrinite  mbe 
loaded  at  the  end  by  the  step  pressure  function  p  =  po  H  (t).  where  po  is  the  ampiimde  of  pressure 
end  K  (t)  is  the  Havyside  function.  Dimersion  effects  in  the  real  configuration  of  the  MDS 
soecimen  and  Kooldnson  tube  have  been  analysed  for  several  idealised  crapezoidai  incident  pulses 
-Sing  a  dynamic  FE  code.  Tnose  caicuiations  revealed  characienstic  features  in  dispersive  changes 
cf  the  transmiaec  pulses  by  the  MDS  specimen. 

•Adiabatic  shear  banding  leads  directly  to  dynamic  fracture  in  Mode  II.  Because  the  stress 
concentrators  are  present  m  every  loading  configuration  this  causes  that  the  dynamic  Mode  11  is  the 
final  Stage  of  the  ASB's.  A  large  study  of  rate  effects  in  Mode  II  fracture  has  been  performed, 
Appendix  n®  1  to  the  Repon. 

Tne  new  experimental  technique  has  been  appiied.  as  an  preliminary  stage,  to  test  an  annea¬ 
led  low  alloy  mild  steel  (O.lTTc  C.  0.589c  Mn).  It  was  attempted  to  rind  conditions  for  catastrophic 
adiabatic  shear  as  a  function  of  imoact  velociry.  up  to  100  irv's.  .An  energy  analysis  have  shown  that 
m  excess  of  certain  impact  velocity,  typically  for  this  steel  Vo  *  90  m'S,  thie  energy  to  break  MDS 
specimen  diminishes  when  the  i.npact  velocity  is  increased.  In  addition,  scanning  micrographs 
were  taken  of  the  fracture  surfaces. 

The  second  senes  of  e.xpcnments  have  been  periorracd  on  quenched  and  temipered  VAR 
■1340  steel.  Again,  the  threshold  of  impact  velocity  has  been  found  for  this  steel  in  excess  of  which 
the  critical  force,  or  energy*  of  fracturing,  diminishes  when  impact  velocity  is  Increased.  The  critical 
impact  velocity  is  estimated  as  Vq  150  m/s.  Scanning  electron  microscopy  revealed  characteristic 
par.ems  on  the  rracture  suriaces  indicating  very  high  temperatures. 
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Extended  abstract 


Since  adiabatic  shear  bands  (ASB’s)  are  frequently  obser\’ed  in  steels  under  high  rate  of 
deformation,  and  they  are  the  dominant  mode  of  failure  in  steels,  the  interest  of  studying  them  is 
sail  actual.  Although  the  ASB's  have  been  previously  studied,  the  effect  of  the  nominal  strain  rate 
as  well  as  the  effects  of  different  impact  velocides  on  ASB’s  development  are  not  throughly  under¬ 
stood. 

A  commonly  used  experimental  technique  in  studying  the  ASB's  is  the  Split  Hopkinson 
Torsion  Bar.  This  technique,  although  quite  effecrive  and  precise,  has  one  disadvantage  that  the 
nominal  rates  of  shearing  are  quite  limited,  typically,  500  s**  <  F  <  2000  S'l,  where  f  is  the  rate  of 
shear  deformation.  In  order  to  expand  shear  testing  to  different  nominal  strain  rates  a  new  experi¬ 
mental  technique  has  been  developed  in  the  Laboratory  of  Physics  and  Mechanics  of  Materials, 
Metz,  France  (LPMM).  .A  modified  specimen  geometry  of  double  shear  has  been  designed  (MDS 
specimen).  The  MDS  specimen  with  flat  faces  is  loaded  directly  by  a  bar  projectile  which  can  be 
accelerated  to  a  desired  impact  velocity  by  a  gas  gun  (1  m/s  <  Vq  <  200  m/s),  where  Vq  is  the 
impact  velocity.  Since  the  specimen  isfcacked  by  a  long  tube  which  transmits  the  force  imposed  on 
the  MDS  specimen  the  shear  stress  history'  can  be  exactly  recorded  by  me  transmitted  elastic  longi¬ 
tudinal  wave.  The  shear  displacement  of  me  MDS  specimen  is  measured  direcdy  by  an  optical 
gage.  Tne  e.xperimencai  scrap  which  consists  of  the  air  gun,  transmitter  tube  (Hopkinson  tube)  and 
optical  displacement  gage  /  e.xtensometer,  permits  for  a  wide  variation  of  the  nominal  strain  rates, 
typically  10^  s'^  ^  F  <  10^  s*^  On  the  other  hand  there  is  no  limitation  of  maximum  shear  strain. 

Since  die  technique  is  new  preliminary  studies,  have  been  undenaken  to  understand  better  the 
mechanics  of  the  test. 

The  MDS  specimen  geometry  has  been  calibrated  with  the  FE.M  and  "ABACUS"  software 
(O.  Oussouaddi). 

Evaluation  of  dispersive  effects  of  elastic  longitudinal  waves  in  rubes  has  been  studied  (J.R. 
Klepaczko  and  S.J.  .Matysiak),  see  Appendix  n’  2  to  the  Final  Technical  Repon.  A  closed- form 
solution  has  been  obtained  for  a  semi-infinite  tube  loaded  at  the  end  by  the  step  pressure  function 
p  =  Po  H  (t),  where  po  is  the  amplitude  of  pressure  and  H  (t)  is  the  Havv’side  function. 

Dispersion  effects  in  the  real  configuration  of  the  MDS  specimen  and  Hopkinson  tube  have 
been  analysed  for  several  idealised  trapezoidal  incident  pulses  using  a  dynamic  FE  code.  Those 
calculations  revealed  characteristic  features  in  dispenive  changes  of  the  transmitted  pulses  by  die 
\fDS  specimen. 

Adiabatic  shear  banding  leads  directly  to  dynamic  fracture  in  Mode  H.  Because  the  stress 
concentrators  are  present  in  every  loading  configuration  and  this  causes  that  the  dynamic  Mode  n  is 
the  final  stage  of  the  ASB's.  A  large  study  of  rate  effects  in  Mode  n  fracture  has  been  performed 
(I.V.  Varfolomeyer  and  J.R.  Klepaczko),  see  Appendix  n®  1  to  the  Final  Technical  Report. 


Different  crack  geometries  in  Mode  II  has  been  analysed,  and  some  of  them  can  be  used  with  the 
new  experimental  setup. 

The  new  experimental  technique  has  been  applied,  as  an  preliminary’  stage,  to  test  annealed 
low  alloy  mild  steel  (0.17%  C,  6.58%  Mn).  It  was  attempted  to  find  conditions  for  catastrophic 
adiabatic  shear  as  a  function  of  impact  velocity,  up  to  1(X)  m/s.  An  energy  analysis  have  shown  that 
in  excess  of  certain  impact  velocity,  or  nominal  strain  rate,  typically  for  this  steel  F  =  1.15  *  10* 
S'‘,  the  energy  to  break  MDS  specimen  diminishes  when  the  impact  velocity  is  increased.  In 
addition,  scanning  micrographs  were  taken  of  the  fracture  surfaces. 

The  second  series  of  experiments  have  been  performed  on  quenched  and  tempered  VAR 
4340  steel.  Again,  the  threshold  strain  rate  has  been  found  for  this  steel  in  excess  of  which  the 
critical  force,  or  energy  of  fracmring,  diminishes  when  impact  velocity  is  increased.  The  threshold 
strain  rate  is  estimated  as  F  =  10^  s*^  Scanning  electron  microscopy  revealed  characteristic  patterns 
on  the  fracture  surfaces  indicating  very’  high  temperatures. 

Although  the  new  e.xperimental  technique  needs  further  improvements  it  appeared  to  be 
useful  and  effective. 


Lis:  of  Keywords :  Steel 

1018  steel 
4340  steel 

High  strength  alloys 
Shear  tests 

.Adiabatic  shear  bands 
Fracture 


-3- 


Introduction 

During  last  decade  a  substantial  progress  has  been  made  in  testing,  analytical  solutions  and 
numerical  calculations  of  the  Adiabatic  Shear  Bands.  It  is  well  known  for  a  long  time,  for  example, 
the  paper  by  H.  Tresca  (1878)  [1],  and  revue  by  W.  Johnson  (1987)  [2],  of  early  works,  that 
metals  can  develop  a  thermal  instability  of  plastic  flow  due  to  heat  generated  during  deformation. 
Zener  and  Hollomon  (1944)  [3],  noted  that  an  increase  strain  rate  is  inevitably  associated  with  a 
change  of  deformation  conditions  from  isothermal  to  adiabatic.  The  same  authors  found  that  the 
plastic  deformation  becomes  unstable  when  the  strain  hardening  rate  (the  tangent  modulus  of  the 
stress-strain  curve)  stans  to  be  negative.  Although  localisation  of  plastic  deformation  in  the  form  of 
shear  bands  is  very  common  in  many  materials,  a  great  deal  of  interest  has  been  recently  devoted  to 
study  the  non-isothermal  cases.  The  present  study  is  also  limited  to  the  non-isothermal  case  when 
the  adiabatic  conditions  of  plastic  deformation  prevail. 

At  moderately  high  and  high  strain  rates  plastic  deformation  in  metals  is  nearly  adiabatic,  as  a 
consequence  the  defomiation  heating  can  lead  to  sizeable  amount  of  flow  softening  and  hence  to 
flow  localization  in  the  lomi  of  ASB,  called  also  the  thermoplastic  catastrophic  shear,  which  can 
intum  lead  to  fracture.  The  ASB’s  are  narrow,  a  small  fraction  of  millimeter,  zones  of  highly  non 
homogeneous  deformation  -ievelopcd  by  a  complicated  interplay  of  strain  hardening,  temperature 
softening,  strain  rate  sensitivity,  dynamic  strain  ageing  and  sometimes  i.ienia  forces.  In  some 
materials,  specially  in  steels  and  titanium  alloys,  a  phase  transformation  or  amorphous-like 
microstructure  may  develop  \vitliin  a  narrow  zone  of  tire  highest  temperatures.  It  may  be  mentioned 
that  there  is  no  agreement  as  to  a  state  of  structure  inside  the  ASB's  for  different  materials.  It  may 
be  mentioned  also  chat  .ASB’s.  with  or  without  phase  transformation,  often  act  as  sites  of  fracture 
initiation  in  Mode  II,  [4].  The  importance  of  the  ASB’s  is  obvious  in  diverse  applications  like 
rolling,  drawing,  machining,  impact  on  structures  and  ballistic  impact. 


.A  number  of  analytical,  numerical  and  experimental  studies  have  been  perfomied  in  attempt 


The  early  revues  on  the  subject  demonstrated  difficulties  in  the  analysis  of  the  problem,  for 
example  A  more  recent  up-to-date  review  was  published  by  Bai  and  Dodd,  [6]. 

Analytical  studies  of  ASB  are  more  numerous  and  they  will  not  be  reviewed  in  this  Rspon. 
Generally,  because  of  simplicity  of  consuiuiiNe  relations,  in  order  to  fmd  a  closed  form  solution, 
inaccepiable  simplislcations  were  automatically  introduced  in  some  ASB  analyses.  More  recent 
anaiy’uc  studies  are  more  acceptable,  for  example  17-101. 

On  the  other  hand,  more  and  nrore  numerical  analy  ses  on  ASB  formation  are  available,  and 
in  general,  the  final  results  are  very  sensitive  to  constitutive  relations  chosen.  Although  very  fre- 
qusntiy  the  fully  rjon-UaRU'  system  of  equations  is  employed  in  such  studies,  she  material  behavior 
characterized  in  the  fonn  of  constitutive  relations  is  quite  simplified.  Such  siraation  may  to 
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Introduction 

During  last  decade  a  substantial  progress  has  been  made  in  testing,  analytical  solutions  and 
numerical  calculations  of  the  Adiabatic  Shear  Bands.  It  is  well  known  for  a  long  time,  for  example, 
the  paper  by  H.  Tresca  (1878)  (1],  and  revue  by  W.  Johnson  (1987)  [2],  of  early  works,  that 
metals  can  develop  a  thermal  instability  of  plastic  flow  due  to  heat  generated  during  deformation. 
Zener  and  Hollomon  (1944)  [3],  noted  that  an  increase  strain  raiC  is  inevitably  associated  with  a 
change  of  deformation  conditions  from  isothermal  to  adiabatic.  The  same  authors  found  that  the 
plastic  deformation  becomes  unstable  when  the  strain  hardening  rate  (the  tangent  modulus  of  the 
stress-strain  curve)  sians  to  be  negative.  Although  localisation  of  plastic  deformation  in  the  form  of 
shear  bands  is  very  common  in  many  materials,  a  great  deal  of  interest  has  been  recently  devoted  to 
study  the  non-isothermal  cases.  The  present  study  is  also  limited  to  the  non-isothermal  case  when 
the  adiabatic  conditions  of  plastic  deformation  prevail. 

At  moderately  high  and  high  strain  rates  plastic  deformation  in  metals  is  nearly  adiabatic,  as  a 
consequence  the  deformation  heating  can  lead  to  sizeable  amount  of  flow  softening  and  hence  to 
flow  localization  in  the  lorm  of  ASB.  called  also  the  thermoplastic  catastrophic  shear,  which  can 
intum  lead  to  fracture.  The  ASB’s  are  narrow,  a  small  fraction  of  millimeter,  zones  of  highly  non 
homogeneous  deformation  ieveloped  by  a  complicated  interplay  of  strain  hardening,  temperature 
softening,  strain  rate  sensitivity,  dynamic  strain  ageing  and  sometimes  inv.ma  forces.  In  some 
materials,  specially  in  steels  and  titanium  alloys,  a  phase  transformation  or  amorphous-like 
microstructure  may  develop  within  a  narrow  zone  of  the  highest  temperatures.  It  may  be  mentioned 
that  there  is  no  agreement  as  to  a  siae  of  structure  inside  the  .ASB’s  for  different  materials.  It  may 
be  mentioned  also  chat  ASB's.  with  or  wit.hout  phase  transformation,  often  act  as  sites  of  fracture 
initiation  in  Mode  I!.  (*i].  The  importance  of  the  ASB’s  is  obvious  in  diverse  applications  like 
roiling,  drawing,  machining,  impact  on  structures  and  ballistic  impact. 

.A  number  of  inalyucal.  numerical  and  experimental  studies  have  been  perfemmed  in  attempt 
to  determine  the  cnocal  conditions  for  the  onset  and  evo'iution  of  castastrophic  thcrmoplasdc  shear. 
The  early  revues  on  the  subject  demonstrated  difficulties  in  the  analysis  of  the  problem,  for 
example  (f  ].  .A  more  recent  up-to-date  review  was  published  by  Bai  and  Doiiki.  [6]. 

.Analytical  studies  of  .ASB  are  nxm  numerous  and  they  will  not  be  reviewed  in  this  Repon. 
Generally,  because  of  simplicity  of  constitutive  relations,  in  order  to  tlnd  a  closed  form  solution, 
inacceptable  simplifications  w-ete  autonuucaily  introduced  in  some  .ASB  analyses.  More  recent 
analytic  studies  are  watt  accepuble>  for  example  (7>  10). 

On  the  other  hand,  more  and  more  numerical  analyses  on  ASB  formation  are  available,  and 
in  general,  the  final  results  are  very  sensttive  to  constitutive  reiatkrns  chosen.  Although  very  fire- 
qu  ?ndy  the  fuUy  non  linear  system  of  equaoons  is  employed  in  such  studies,  die  material  bduvior 
characterized  in  the  focm  consittutive  relations  is  quite  simplified.  Such  siraadon  may  lead  to 
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some  misinterpretaiions  of  how  ASB’s  evolve  as  a  function  of  both  the  initial  and  boundary 
conditions. 

The  state  of  the  an  is  that  relatively  large  quantity  of  theoretical  studies  have  been  so  far 
published  with  much  less  effon  put  into  experiment.  Most  of  experimental  studies  make  use  of 
Split  Hopkinson  Torsion  Bar  (torsional  Kolsky  apparatus),  [11-13].  They  are  based  on  fast  torsion 
of  thin  tubular  specimens  of  shon  lengths.  Such  test  with  SHTB  has  a  very  limited  range  (in  the 
logarithmic  scale)  of  the  nominal  shear  strain  rates,  usually  around  10^  s'\  and  the  boundary 
condidons  imposed  during  the  test,  although  well  defined,  are  not  constant,  [14].  Thus,  majority 
of  experimental  studies  were  Limited  to  a  thin  tubular  geometry'  and  to  the  nominal  strain  rates 
around  lO-’  S‘^  The  other  range  of  rates,  albeit  much  higher,  is  met  in  external  ballisdcs  and 
explosive  loading.  .Although  it  is  very’  easy  to  produce  a  netxvork  of  ASB’s  during  ballisdc  or 
explosive  loadings,  they  do  no:  provide  a  good  e.xperimental  basis  to  srudy  fundamentals. 

The  most  imponant.  and  so  far  unresolved  class  of  problems,  is  the  effect  of  inidai  and 
boundarj.-  conditions  on  the  onset  and  formation  of  the  .ASB’s.  From  the  point  of  view  of  experi¬ 
ments  the  followarg  loading  schemes  cn  be  specified : 
i.  inidaily  uniform  defomaaon  field  with  small  inical  penurbadons ; 
the  Lnitial  penurbadons  w  hich  are  usually  issumed ; 

a.  small  uefict  of  geomeuy ; 

b.  small  uiemiai  pemurbations  i 
c  ’neat  sinks. 


Imecscd  deformation  fries  w  i:h  unsantaneous  stress  or  smiin  co.ncenmito.'s ; 


a.  conc-o*ied  Icacing  conenons ; 


b.  dmec:  umpac:  a:d:::eren:  veiocises  tinciuong  ballisdc  umpacti  and  dutereat  projeeme; 


target  geemem 


$ 


:;i.  E.xpiO»;xe  loading ; 

a.  e.\par,s;on  o:  rubes  ifragme.nution  .ASB’si ; 

b.  shock  waves  generated  or.  limited  sunaces. 
tv  Conmoiled  shock  wav'cs ; 

a.  piate.-piaw  impact  ^micro  .AS 8  si. 

V .  Meal  foctning  processes ; 

a.  machir.tne ; 

b.  high  spee  d  reeling . 

c.  high  speed  dnwtng. 

The  systemadcs  shown  ^ve  points  out  on  a  wtde  class  of  loading  condidons,  and  in  turn,  on  a 
wide  class  of  inioal  and  boundary  conditions  under  which  the  ASB's  can  be  generated.  The  cases 
iii..  iv.  and  v.  will  not  be  discussed  here,  however,  the  most  fundamental  cases.  L  and  it.  are 
wonh  of  discussion.  The  case  i.  is  clcariy  relued  to  SHTB  techniqiK.  U  can  be  point  out  that  the 
aenve  specimen  length  and  minai  condinons  can  inf,ue.^ct  the  outcome  in  the  fwm  of  .ASB.  For 


i 
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example  Litonski  [15]  has  furnished  a  detailed  analysis  of  the  deformation  in  torsion  of  a  thin- 
walled  tube  with  an  initial  geometric  imperfection  in  the  form  of  a  slightly  thinned  section.  Tlie 
analysis  was  repeated  in  a  similar  way  by  Costin  et  al  [11]  but  with  comparison  to  SHTB 
experiments  on  1018  CRS  (cold  rolled  steel)  and  1020  HRS  (hot  rolled  steel).  Sometime  later 
Litonski  [16]  have  shown  by  a  similar  numerical  analysis,  but  the  thin- walled  tube  without 
geometric  imperfection  has  been  analysed  and  with  the  heat  sinks  at  the  ends,  that  the  ASB  occurs 
above  cenain  critical  strain  rate  of  the  order  1.0  s*L  Another  numerical  study  was  reponed  in  [17] 
where  a  thin-walled  tube  of  a  constant  cross-section  had  a  Gaussian  imperfec-tion  of  the  yield 
stress  Xy.  i-C-  'Ey  (x)  =  tyo  G  (x),  where  G(x)  is  the  Gaussian  distribution  and  tyo  is  the  mean  yield 
stress.  It  was  shown  in  this  study  that  at  lower  nominal  strain  rates,  of  the  order  0.1  s'^  the  strain 
and  temperature  fields  differ  for  the  adiabatic  and  heat  conduction  cases.  When  the  nominal  strain 
rate  was  increased  to  1 S'^  the  differences  were  substantially  reduced. 

Also  later  the  geometry  imperfections  were  the  most  frequently  studied  cases.  In  [18]  are 
reported  numerical  simulations  of  ASB’s  observed  on  1020  in  [11].  Plastic  instability  and  flow 
localization  in  shear  at  high  rates  of  deformation  were  studied  in  [19]  by  assuming  a  geometrical 
imperfection  in  specimen.  In  [20]  the  effects  of  material  imperfections  via  geometric  imperfections 
on  flow  localization  in  SHTB  test  were  studied. 

Shawki  and  Clifton.  [21]  presented  a  number  of  analytic  solutions  and  numerical  studies, 
including  torsion  of  a  thin  tube  with  geometrical  imperfections,  with  one  or  more  circumferential 
grooves  and  with  themtally  isolated  sides.  Three  simple  constitutive  relations  were  used  in  those 
calculations.  It  was  found,  as  expected,  that  at  high  .nominal  strain  rates  of  the  order  10^  to  10^  s'^ 
the  localization  mechanism  which  accounts  for  .ASB  formation  is  exclusively  adiabatic.  For  this 
mechanism,  at  adv.inced  deformation  the  thermal  softening  dominates  in  reduction  of  stress  and  an 
increase  of  local  shear  strain  leading  ic  catastrophic  local  shearing.  .At  small  strains  a  positive  strain 
hardening  dominates  and  localization  does  not  occur.  The  rate  of  localization  depends  strongly  on 
the  strain  rate  sensitivity.  The  qualitative  discussion  as  offered  above  is  consistent  with  earlier 
findings  and  physical  intuition.  The  qualitative  analyses  of  the  ASB’s  formation  for  particutar 
metals  and  alloys  and  for  different  boundary  conditions  are  sull  waiting  to  be  solved. 

.A  more  recent  numerical  study.  (22)  base  shown  that  the  heat  sinks  introduced  at  the  ends  of 
a  thin- walled  tube  can  modify  the  whole  process  of  the  ASB  formadon  in  comparison  to  the 
assumption  of  the  adiabatic  boundary  conditions.  In  addition  to  the  heat  sinks  a  temperature 
perturbation  was  introduced.  .A  small  local  maximum  in  the  initial  temperature  distribution  along 
the  tube  axis  developed  fonnatioo  of  the  ASB  just  in  this  cross  section.  A  simple  linear  Anhenius 
model  combined  with  a  power  strain  hardening  was  employed  in  those  calculations,  the  nominal 
strain  rate  was  assumed  as  1.6  x  1(P  s'L  Again,  numerical  results  discussed  above  clearly  indicate 
the  importance  of  the  initial  and  boundary  conditions  on  flnal  geon^try  and  formaiioo  of  the 
ASBs. 
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Transition  from  isothermal  to  adiabatic  deformation 

The  localizadon  of  deformation  in  the  form  of  ASB’s  is  observed  in  thin-walled  tubes  when 
the  nominal  strain  rate  exceeds  cenain  value.  Obviously,  a  transition  exists  between  pure 
isothermal  and  pure  adiabatic  regimes  of  deformation.  This  transition  depends  in  the  first  place  on 
geometry  of  deformed  body  and  on  effectiveness  of  heat  evacuation  from  the  heat  zones.  When 
correct  boundary  and  initial  problems  are  posed  it  is  possible  to  estimate  the  range  of  strain  rates 
within  which  the  transition  occurs.  Preliminary  calculations  to  find  an  approximate  strain  rate  range 
of  transition  were  reported  in  [15]  and  [17].  In  the  second  case  numerical  calculations  were 
performed  for  the  length  of  tube  L  =  30  mm  and  initial  temperature  673K,  and  three  cases  were 
studied  :  isothermal,  adiabatic  and  mixed,  when  it  was  allowed  for  heat  conduction.  Three  nominal 
strain  rates  were  assumed  :  0.01  ;  0.1  and  1.0  S‘h  The  final  results  of  those  calculations  have 
shown  that  the  difference  between  adiabatic  and  mixed  cases  is  substantially  reduced  for  shear 
strain  rate  1.0  s'^  A  similar  calculation  was  performed  in  [16]  for  a  thin-walled  tube,  L  =  5.0  mm, 
and  with  the  heat  sinks  at  the  ends,  again  the  initial  temperamre  was  assumed  To  =  673K,  and  four 
values  of  nominal  strain  rates  were  considered  :  0.01  ;  0.1  ;  1.0  and  10.0  S'^  Similar  plots  were 
produced  for  the  complete  case  (heat  conduction  included).  It  was  found  that  the  transition  must 
occur  between  nominal  strain  rates  0.1  and  1  s*^ 

A  more  exact  transition  analysis  for  thin-walled  tubes  of  different  lengths  was  performed 
more  recently  by  Oussouaddi  and  Klepaczko  [23].  Tne  finite  difference  technique  has  been  applied 
with  relatively  exact  constitutive  relation,  and  calculations  were  carried  out  at  different  nominal 
strain  rates  in  shear  from  1.0  s**  to  10^  s*^  at  To  =  300K.  In  addition,  the  effect  of  length  on  the 
isothermal/  adiabatic  transition  was  studied.  Since  the  heat  sinks  were  assumed  at  the  tube  ends  the 
transition  was  denned  as  the  maximum  of  partial  derivative  of  the  temperature  gradient  with  respect 
to  the  nominal  strain  rate 

s  _  3  /5T] 

SlogfJ^I  (1) 

where  x  is  the  axial  coordinate  ;  -  L/2  <  x  <  L/2,  T  is  the  absolute  temperature  and  Fn  =  dTn/dt  is 
the  nominal  strain  rate  in  shear,  Fn  =  (r/L)  tp,  r  is  the  mean  radius  and  tp  the  angle  of  twist.  The 
temperature  gradients  as  a  function  of  log  Fn  for  L  =  2.0  mm  are  shown  in  Fig.  la,  whereas  Fig. 
lb  shows  variation  of  c  also  as  a  function  of  log  F-.  It  is  clear  that  the  transition  can  be  well 
defined  by  the  maximum  value  of  q.  The  critical  values  of  (Fnlc-  are  given  in  Table  1,  and  also 
shown  in  Fig.  2  as  a  function  of  L.  The  results  suppon  qualitatively  the  physical  intuition  that  the 
critical  strain  rate  increases  in  proportion  to  the  thermal  conductivity,  which  is  the  lowest  for  steel 
and  the  highest  for  copper.  When  L  decreases  from  10  mm  to  1.0  mm  the  critical  strain  rate 
increases  approximately  twofold  for  all  three  material  which  were  smdied. 
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Those  results  clearly  indicate  the  imponance  of  the  boundary  and  initial  conditions,  in  both 
experiment  and  numerics,  not  only  in  the  case  of  transition  isothermal/adiabatic,  but  also  in  for¬ 
mation  of  ASB’s. 

Table  1 

Critical  shear  strain  rate  (s'O 


Lenssth  (mm) 

Steel 

Aluminium 

CoDPer 

1 

58.0 

89.0 

112.0 

47.5 

68.5 

84.5 

4 

40.5 

54.0 

66.0 

6 

36.0 

50,0 

59.0 

10 

31.0 

42.0 

49.0 

Shear  strain  of  instability 


The  earliest  studies  on  the  onset  of  the  shear  localization  were  limited,  by  analogy  to 
Considere  condition  of  instability  in  tension  test,  to  the  load  instabilit>'  dM/d(p  =  0.  The  condition 
dM/d(p  =  0,  where  M  is  the  torque  and  o  is  the  angle  of  twist,  gives  possibility  to  determine  the 
critical  instability  strain  if  a  specific  constitutive  relation  is  assumed.  It  was  confirmed  by  Litonski 
[15],  by  numerical  application  of  the  flow  localization  model  (Marciniak-Kuczynsld  approach 
[24]),  that  the  onset  of  instability  corresponds  to  the  maximum  load,  that  is  the  condition  dlvl/d9  = 
0  was  confirmed  Cj  ihs  onset  of  instability.  Many  authors  derived  formulas  for  the  critical  strain 
using  different  empirical  consututive  relations,  a  review  of  those  derivations  was  published  in  [25]. 
Experimental  confirmation  of  the  onset  condition  dM/dtp  =  0  was  provided  for  a  low-alloy 


structural  steel  in  [13].  .A  high  speed  photography  combined  with  SHTB  technique  was  applied  in 
those  study. 

Since  up  to  the  maximum  of  torque  the  shear  deformation  is  uniform  the  condition  for  load 
instability  reduces  to  the  formula 
dx 


•  =  0 


dT 


(2) 


Because  many  authors  use  constitutive  relations  in  the  form  of  multiplication  function 


t=fiCDf2(r)f3(n  (3) 

a  more  general  analysis  of  the  condition  (2)  will  be  given,  KIcpaczko,  [26].  If  the  mechanical 
equation  of  state  is  assumed  and  the  history  effects  arc  neglected  one  can  write 


then  the  condition  (2)  leads  to 


.orjir"  'oTlnf  ^  '  [^rjTr  dT  (5) 

Condition  (5)  can  be  satisfied  only  in  very  specific  conditions  of  deformation.  One  of  them  is  the 
adiabatic  process  of  deformation,  then 

£L  = 

dT  dr|ADLJ£.\T!C  (6) 

Condition  (S'!  can  be  analysed  for  variety  of  deformation  histories  f  (F)  and  T  (F),  but  the  sim- 
ciest  case  is  usually  iiimted  to  tne  constant  strain  rate,  in  general 

^1  1^  '  (  il  !  =n 

l^flrldriA  (7) 

It  is  important,  then,  to  have  an  estimate  of  the  adiabatic  increase  of  temperature  T.a,  (F)  and 
(dT/dF)A  due  to  plastic  work  convened  into  heat.  Assuming  constitutive  relation  (2)  it  follows 

from  the  principle  of  energ>'  conser.’ation 

« 

cT  [I  •  C  CTF)]  f  1 CT)  f:  (Q  ^3  (^) 

dF  A  pCTlCpCD  (8) 

where  (T.F)  is  a  ccefficient  taking  into  account  the  stored  energy  in  the  matetial,  p  (T)  is  the  mass 
density  of  the  material  and  C;  (T)  is  the  specific  heat  at  constant  pressure,  both  p  and  Cp  are 

function  of  temperature. 

The  partial  differennals  can  be  found  after  (3)  as  follows 


Introduction  (9)  into  condition  (4)  yields 


After  elintinaticn  of  (dT/dF)A  using  (8)  the  condition  (10)  is  transformed  into  the  following  form 
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faf. 


af, 


filrj  lar, 


+  f2(rjf3(f)  - 


1-^ 


/afi\ 


ar;  p(T)CpCD  f 


'dr 


ar;  Idr 


=  0 


where  H  is  the  strain  at  the  onset  of  instability. 

If  the  strain  rate  is  constant,  (df/dF)  =  0,  condition  (11)  reduces  to 


^2 


i-C 


■  =  0 


p(T)Cp{T) 

Since  only  f:  (Fc)  depends  on  the  critical  strain  the  explicit  equation  for  fj  (Fc)  is 

■  -.1/2 

(af^/ar)  pCDCpCia 


f2(r.)= 


f3(F)(l-0 


01) 


(12) 


(13) 


The  expression  has  real  and  imaginary  part.  Inversion  of  f2  (Fc)  makes  it  possible  to  find 


rc=f2 


A(r,r,T 


1/2 


(14) 


where  -  A  is  the  expression  in  the  square  brackets  of  (13).  Existence  of  the  real  f2  (Fc)  is  possible 
only  if  A  is  negative.  Since  p,  Cp,  (1-^)  must  be  always  positive  and  the  function  of  strain  rate 
sensitivity,  U  (O.  is  assumed  in  this  case  also  positive,  the  only  term  which  may  be  negative  is 
(afa/aF)  /  (afi/ai).  The  most  common  case  is  the  thermal  softening  which  leads  to  negative  value 
of  afi/aT,  of  course,  if  dtVaF  is  at  the  same  time  positive.  Another  possibility  is  if  the  tangent 
modulus  afyar  is  negative  and  at  the  same  time  there  is  no  thermal  softening.  The  role  of  a 
positive  rate  sensitivity  is  quite  interesting,  that  is  if  t's  (O  is  an  increasing  function  of  strain  rate  F. 
The  positive  rate  sensitivity  has  a  negative  effect  on  the  onset  of  adiabatic  instability,  that  is  Fc  is 
reduced  when  strain  rate  is  increased.  The  positive  rate  sensitivity  increases  production  of  plastic' 
work  convened  into  heat.  However,  after  the  critical  strain  Fc  is  reached,  F  >  Fc,  a  positive  rate 
sensitivity  diminishes  local  strain  gradients  in  ASB’s.  Since  the  assumption  was  made  (dF/dF)  = 
0,  function  f}  (H  plays  a  role  of  parameter,  it  must  be  constant  but  rate  dependent. 

In  many  publications  constitutive  relations  are  employed  which  fail  into  the  class  of  equation 
(3).  For  example,  in  [15-17]  and  other  publications,  the  following  relation  was  used  in  numerical 
stU'iies 


x  =  B(l-aT)(l 


(15) 


where  n  is  the  strain  hardening  index  and  m  is  the  rate  sensitivity,  a,  b,  B  are  empirical  constants. 
In  this  case  the  thermal  softening  term  has  been  linearized.  Another  authors  employed  the  follo¬ 
wing  constitutive  relation,  for  example  [7, 9, 21) 
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t  =  B 


r^o 

It 


F^ 

0 

m 


(16) 


where  B,  v,  n,  m  are  respectively  plasticity  modulus,  temperature  index,  strain  hardening  index 
and  logarithmic  strain  rate  sensitivity,  Tq,  Fq  and  Fq  are  normalization  constants.  Relation  (16) 
falls  in  the  class  of  non-linear  liquids  since  r  =  0  if  F  ->  0.  This  relation  was  used  many  times  for 
both,  numerical  solutions  and  penurbation  analyses,  for  example  [7,  21,  27,  28],  and  others. 

To  discuss  funher  the  condition  of  adiabatic  instability  it  is  assumed  tiiat  the  flow  stress  is 
represented  by  a  more  specific  relation 


n  .  m 

z=B(r)r  r 


(17) 


where  B  (T)  is  known  temperature-dependent  plasticity  modulus.  The  condition  (10)  yields  the 
following  relation 


n-i  ,  r 

n  r„  -  F  r 


j  o3  j 
loT.' 


1-C 


(18) 

it  is  assumed  that  the  mass  de.nrity  po  and  the  coefficient  of  stored  energy  remain  constant.  After 
(IS)  the  critical  shear  strain  at  t.ne  onset  of  instabilitv’  is 


tp  C:,vTi 


ri-OF 

'  Cl  ' 


r>-l 


(19) 


•A  positive  role  of  strain  hardening  in  delaying  the  adiabatic  instability  was  well  documented  in 
early  oublications.  Since  tne  specific  heat  is  an  increasing  function  of  temperature  it  also  delays 
instability.  The  most  impo.-tar.;  factor  is  the  mathematical  form  of  the  (-  5B/cT),  especially  in  the 
case  of  post-crincal  analyses. 

.Adiabatic  increments  of  temperatu.’s  at  the  point  of  instability,  9”'0r  =  0.  are  not  usually  very  high, 
and  constancy  of  n  and  m  are  acceptable.  Constant  values  of  n  and  m  are  only  the  first 
approximation  for  sm.al!  increments  of  temperature.  Since  it  has  been  shown  in  [29]  and  later  in 
[  1?)  that  the  rise  of  the  post-cr.ncal  temperature  for  steels  is  of  the  order  425  K  S  AT  S  595  K  at 
the  centra!  pan  of  tiie  ASB,  the  constancy  of  n  and  m  is  a  ver'  crude  approximation.  .A  similar 
conclusion  was  reached  in  [301  by  numerical  analysis  of  a  thin-walled  torsion  test  spcciman.  The 
thermal  softening  -  (dB/oT)  was  assumed  linear  in  the  first  calculations  of  the  post-critical 
deformation  of  copper  and  bilinear  in  the  second  one.  Substantial  differences  were  revealed  in  the 
evolution  of  stresses  in  both  stages  of  deformation,  up  to  the  instability  points  and  during 
localization.  Another  paradox  of  the  temperature-independent  Icgarithmic  rate  sensitivity,  if  the 
thermal  softening  is  negative,  lies  in  the  fact  that  the  stress  differences  at  two  consunt 

strain  rates  increase  according  :o  formula  ( 17)  when  tcmperatu.'t  is  decreased,  this  is  in  complete 
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contradiction  to  the  real  behavior  of  metals,  [31].  A  more  rational  approach  is  to  use  some  form  of 
Arrhenius  relation  to  take  into  account  at  the  same  time  thermal  softening  and  temperature- 
dependent  rate  sensitivity,  for  example  [21]  and  also  [22]. 

Because  of  complexity  in  formation  of  the  ASB’s  the  analytic  studies  in  this  area,  although 
imponant,  have  qualitative  character,  mostly  due  to  simplistic  constitutive  relations.  The 
importance  must  be  focused  on  a  full  numerical  approach  with  physically-based  constitutive 
relations. 


Some  results  of  numerical  studies 


Among  publications  on  numerical  analyses  of  ASB’s  only  few  take  into  account  the  fully 
non-linear  system  of  governing  equations.  However,  the  crucial  problem,  even  for  the  complete 
formulation,  is  die  choice  of  constitutive  relations. 

Two  boundary  value  problems  are  the  most  interesting  from  the  point  of  view  of  experi¬ 
mental  techniques : 

i.  fast  torsion  of  a  thin- walled  tube  (interesting  due  to  SHTB  technique) ; 

ii.  fast  shearing  of  a  layer  of  finite  hight  (interesting  due  to  a  new  experimental  technique  descri¬ 
bed  in  the  next  pan  of  this  repon). 

Relatively  complete  analysis  of  torsion  was  reponed  in  [22].  Numerical  smdy  of  a  layer  as  a 
complete  non-linear  problem,  fully  coupled  w’ith  temperature,  was  reported  in  [32].  The  fully 
coupled  problem  was  formulated  including  heat  conduction,  changes  of  specific  head  as  a  function 
of  temperature,  complete  effects  of  temperature  on  mechanical  properties  and  inertia.  The  layer 
with  geometric  imperfection  was  assumed,  and  the  top  and  bottom  of  the  layer  were  adiabatically 
isolated.  A  complete  phenomenological  constitutive  relation  has  been  used  in  those  calculations.  A 
more  exact  analysis  and  discussion  of  this  constitutive  relation  is  published  elsewhere,  [33].  The 
fully  temperature  coupled  constirutivc  relations  are  as  follows 


t=B(e)(r.-r)  r  +<  ii(e)|r-rj  >  po, 

where  6  is  the  homologous  temperamre,  0  =  T/Tm.  and  Tm  is  the  melting  temperature  of  a  material 

•  • 

under  consideration,  q  is  the  pseudo-viscosity  and  <  •  >  is  the  operator,  <->  =  0ifr^ro  and 
<  •  >  =  I  if  r  >  To,  To  is  the  threshold  strain  rate.  Both,  strain  hardening  index  n  and  logarithmic 
rate  sensitivity  m  are  coupled  with  temperature.  For  BFC  metals  u  and  n  aie  expressed  as  a 
function  of  6  as  follows 


o(e)  =  no— 


(21) 
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m  (0)  =  a  6  exp  (-  b  6)  4-  6  ^27) 

where  rio  is  the  strain  hardening  index  at  OK,  and  a,  b,  tto  and  r  aie  constants,  )i  (0)  is  tempera¬ 
ture-dependent  shear  modulus  of  elasticity 


■* 

1 

1  -  0  exp  j 

1 

Si 

i  elij 

(23) 


where  po  is  the  shear  modulus  at  OK  and  6u  is  a  constant.  All  constants  for  two  steels  (1018  CRS 
and  1020  HRS)  are  identified  in  the  paper  which  was  mendoned  above.  Variation  of  specific  heat 
was  taken  into  account  according  to  Debye  formulation.  Geometric  imperfection  was  assumed  in 
the  form 


V'  (y)  =  w 


1  +5sin 


(24) 


where  Wq  is  the  layer  thickness  and  h  i-*  the  layer  height,  5  and  e  are  geometry  parameters.  Equa¬ 
tions  of  momentum,  b'  lance  o:  energy  with  adequate  Fourier  constant  for  heat  conduction  and 
compatibility  condition  were  assumed  iii  their  standard  form  in  those  calculations.  Numerical 
calculations  were  performed  using  the  finite  difference  implicit  method.  The  duration  of  the 
successive  time  increments  were  chosen  ir  such  a  way  that  the  algorithm  was  unconditionally 
stable.  The  boundaty  ccntiicions  were  assuiued  that  the  bottom  sunace  of  the  layer  was  fixed  and  a 


constant  velocity  was  imposed  to  the  top,  a  steady  smce  process. 

The  input  for  those  calculations  was  the  nominal  strain  rate  in  shear  Tn  =  V/h.  The  output  of 
calculations  consisted  o.'  a  detailed  spati‘’J  history  or  all  important  phys*cai  quantities  like  shear 
stress  T  (y,t),  strain  f  (y,t),  strain  rate  f  (y.t).  absolute  temperature  T  (y.t'i  ai.i  local  velocity  v 
tv.ti.  Thus,  the  simulation  provided  a  wnole  dynarrau  nistory  of  the  A3B  development  at  different 
nominal  strain  rates  Tn- 

Two  values  of  critical  nomirai  strains  were  evaluated,  the  instability'  strain  fc  associated  with 
the  mxvimum  force  and  the  final  localization  strain  Tyx  defined  as  follows 


''n-'^LOC  limr5  =  »c 

• 

(25) 

v/here  Pa  is  the  local  sti-ain  rate  ir.  the  strongest  section  of  the  layer  and  Fb  is  the  local  strain  rate  in 
ire  cross  section  v  here  localization  occurs.  N’umcricai  solutions  obtained  for  two  steels  within  the 
range  of  nominal  strain  rates  10^  s*'  S  Fn  SI  10*  S’*,  (1018  CRS  and  1020  HRS)  have  shown  a 
minimum  of  both  values  of  Fn,  i.e  Fc  and  Floc  Oie  nominal  strain  rate  -  2x10^  S''.  Values  of 
Flxx:  were  obtained  numerically  using  the  condition  Floc  =  lO-  Fn-  In  the  range  of  Fn  lower  than 
Fn  =  2x10^  S"^  both  nominal  strains  F®  and  Flcx:  decrease  when  fn  is  increased,  the  trend  in 
accor(ia.ice  to  formulas  (13)  and  (19).  In  the  range  of  the  nomin?l  strain  rates  higher  than  Fn  * 
2x10^  s*‘  a  slight  increase  of  F.-  and  Floc  is  observed.  Since  calculations  were  performed  twice 


without  pseudo-viscosity,  T]  =  0,  and  with  pseudo-viscosity  the  positive  effect  of  excessive  strain 
rate  sensitivity  on  the  localization  is  obvious.  The  ratio  Floc/Tc  increases  at  the  higher  rate  if  the 
pseudo- viscosity  is  included  into  calculations,  if  only  the  logarithmic  rate  sensitivity  is  taken  into 
calculations  the  ratio  Floc/Fc  increases  only  slightly  as  a  function  of  Fn.  It  has  been,  thus,  shown 
in  the  paper  discussed  above  that  for  a  given  material  and  a  given  geometric  imperfection,  there 
exists  a  domain  of  nominal  strain  rates  Fn  where  development  of  the  catastrophic  shear  in  the  form 
of  the  ASB  is  the  most  likely  to  occur.  A  similar  minimum  of  Fc  was  also  found  in  [34]  using  a 
finite  element  technique,  and  confirmed  later  in  [35]  by  an  approximate  method.  In  both  cases  the 
initial  temperature  defect  was  assumed.  Attempts  to  find  the  final  thickness  of  ASB’s  after 
numerical  results  reponed  in  [32]  were  not  quite  successful.  Although  it  was  found  that  ASB’s 
were  formed  the  asymptotic  width  was  difficult  to  estimate  since  at  Floc  =  10^  Fn  the  evolution  of 
ASB’s  were  still  not  finished.  One  of  reasons  why  it  was  difficult  to  specify  a  finite  thickness  of 
the  ASB  is  the  excessive  strain  hardening  introduced  by  the  term  at  large  strains,  even  when 
the  correct  relation  for  n(T)  was  used. 

As  an  result  a  completely  new  series  of  numerical  calculations  were  undertaken,  but  Lhis  time 
a  physically  based  set  of  constitutive  relations  have  been  used.  Here  only  some  results  will  be 
discussed,  a  more  thorough  of  that  study  is  given  elsewhere,  [36].  The  constitutive  modeling  is 
represented  by  a  consistent  approach  to  the  kinetics  of  macroscopic  plastic  behavior  of  metals  with 
BCC  and  FCC  structures,  [37, 38].  The  constitutive  formalism  is  used  here  with  one  state  variable 
which  is  the  total  dislocation  density  p.  The  notion  is  adopted  that  plastic  defonnation  in  shear  is 
the  fundamental  mode  of  deformation.  It  is  assumed  that  at  constant  microstructurc  the  flow  stress 
T  has  two  components :  the  internal  stress  and  the  effective  stress  x*,  thus 


x  =  x., 


h(f.Tl.r.TsTT^-rX*(r.T)s 


‘JSTR  ^  ‘  •  ^.'STR  (26) 

where  h  (f,T)  arc  the  histories  of  plastic  deformation,  F  and  f  arc  respectively  plastic  strain  and 
strain  rate.  The  internal  stress  x^  is  developed  by  long  range  strong  obstacles  to  dislocation  mo¬ 
tion,  and  the  effective  stress  x*  is  due  to  thermally  activated  short  range  obstacles.  In  fact,  the 
internal  stress  x^  must  be  also  rate-  and  temperature-dependent  via  dynamic  recovery  or 
annihilation  processes.  To  describe  completely  strain  rate  and  temperature  effects,  including 
evolution  of  microstructure  and  smtin  hardening,  an  evolution  equation  must  be  defined.  In  this 
case  a  simple  evolution  equation  has  been  adopted  which  is  adequate  to  account  for  strain 
hardening  in  mild  steels.  [37] 


dT  ^  (27) 

whem  po  is  the  initial  dislocation  density  and  k.  (T.T)  is  the  annihilation  factor.  The  multiplication 
facnx'  Mq  is  related  to  the  mean  free  path  X  of  dislocation  storage  by  the  fonnula  Mn  ■  1/bX,  where 
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b  is  the  modulus  of  Burgers  vector.  For  constant  temperature  and  strain  rate  the  explicit  form  for  p 
can  be  found  after  integration  of  (27) 


If  strain  rate  and  temperature  are  not  constant,  like  in  the  case  of  ASB  formation,  the  evoludon 
equariot:  must  be  numerically  integrated  to  find  current  value  of  p.  If  p  is  found  the  internal  stress 
is  also  found  from  the  reladon 
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Pm  =  fp 


(34) 


where  Go  is  the  energy  of  the  double  kink  formation  and  x  is  the  Peierls  stress  at  OK,  a*  is  the 


activation  distance,  a*  =  n  b,  [37, 41].  The  frequency  factor  Vk  is  assumed  in  the  form 


where  fo  is  the  initial  fraction,  eq.  (34),  pm  is  the  mobile  dislocation  density,  Vd  is  Debye  fre¬ 
quency,  fj,  T,  are  constants.  The  set  of  equations  from  (32)  to  (35)  provide  complete  evolution  of 
the  effective  stress  x*. 

This  constitutive  modeling  involves  two  kinds  of  constants,  the  absolute  physical  co.istants 
and  the  constants  also  based  on  physics  but  specified  for  a  particular  metal  or  alloy.  The  complete 
discussion  of  all  constants  will  be  given  elsewhere,  [36],  but  some  basic  values  are  given  below  : 
Mn  =  1.2x1010  cm  - ;  T.  =  558  K  ;  Tr  =  300  K  ;  iHo  =  1.47x10^,  To  =  10*  s  ' ;  Go  =  0.315  c  V  ; 
ko  =  9.23  ;  fi  =  7  ;  =  520  .VIPa  ;  n=2  ;  fo  =  0.05  ;  po  =  6.2x108  cm  - ;  q  =  4/3  ;  p  =  3/4.  Those 

constants  are  represenianve  for  a  mild  steel  XC18  (French  Standards)  after  annealing,  0.17%  C ; 
0.58%  Mn;  0.21%  Si. 

Numerical  calculations  were  performed  using  the  same  finite  difference  procedure  depicted 
previcusly,  but  the  hight  of  the  layer  was  assumed  h  »  2.0  nun  and  the  number  of  elementary 
sublayers  was  assumed  100,  so  that  the  elementary  sublayer  had  the  thickness  Ah  »  20  lun.  With 
such  value  of  Ah  an  analysis  of  the  terminal  thickness  of  the  ASB  could  be  attempted. 

Results  of  seventeen  computer  simulations  for  different  nominal  strain  rates,  80  S'l  S  1 
4.5x10*  s'l,  will  be  bnefly  presented  here.  The  nominal  instability  strain  Pc  and  the  nominal 
localization  strain  Ftoc  are  shown  as  a  function  of  log  r„  in  Fig.  3.  The  minimum  of  the  instabi¬ 
lity  strain  Pc  is  found  again  in  the  region  of  nominal  strain  rates  above  10^  s  ^  log  P^  ■  2.3,  Pq  ■ 
200  S'^  This  is  better  visible  in  Fig.  4  with  expanded  scale  of  Pc,  The  stabilizing  effect  of  strain 
rate  and  perhaps  inertia  are  cteaiiy  shown  in  Fig.  S  where  the  ratio  P^loc  i-^  plotted  as  a  functitm 
of  log  Pc-  The  ratio  begms  with  value  around  2.0  at  »  10^  s**  and  ends  with  -  6.0  at  Pn  ■  10* 
s^  Around  the  range  of  the  nominal  strain  rates  *>5x10*  the  maxiinum  temperature  at  the  center 
of  the  ASB  is  found  to  be  quite  close  to  the  melting  point  for  iron,  Tm  »  17SS  K.  This  is  why  the 
last  point  of  localization  strain  Ptoc  calculated  for  the  highest  Pg  is  (ki  the  same  level  in  Fg.  3  as 
the  preceding  point,  Ptjoc  ■  2.5.  A  high  probability  that  at  the  latest  stages  of  ASB  fbrmatioo  the 
lemperaniie  at  the  central  pans  ASB's  exceeds  the  melting  point  was  suggested  in  [42].  It  was 
found  after  examination  of  fracture  surfaces  formed  during  reosUe  lesu  of  hi^  strength  steels  that 
the  tops  of  the  dimples  separated  by  a  fast  shearing  are  spheroidized,  appaitaily  from  melting 
temperature. 

The  most  sensitive  paraioeier  for  numerical  analyses  of  the  ASB's  is  die  final  width  of  the 
band.  In  the  present  analysis  a  gradient  method  was  assumed  to  defuse  the  final  width  of  shear 
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band.  The  profiles  of  shear  straLn  F  (y)  and  temperature  T  (y)  as  shown  in  Fig.  6  and  Fia.  7  for  Fn 
=  10-'  s-  and  for  differen:  ieveis  of  the  nominal  strains  F„,  ser\'ed  as  a  basis  for  all  strain  rates 
anaivsed  to  nnd  strain  and  temperature  gradients  gr  (y)  and  gj  (y)  at  different  levels  of  the  nominal 
s~ain.  The  strain  gradients  gr  iy)  obtained  after  such  procedure  for  Fn  =  10^  s*^  at  different  levels 
of  nominal  strain  are  shown  in  Fig.  8.  The  width  of  the  ASB  has  been  defined  as  an  asymptotic 
value  of  the  distance  ber.veen  the  gradient  maxima, 

lim  Ay_  =  2a  or  iim  Ay^  =  2a  ^^6) 

F  — ^  Floc  T  — *  Tloc 

where  Aym.  is  the  distance  benveen  the  gradient  maxima,  both  for  strain  and  temperature.  Fig.  9 
shows  details  how  the  asymptonc  half-value  a  of  the  ASB  thickness  was  determined.  The  dashed 
line  shows  terminal  evolution  of  the  band  thickness.  The  same  procedure  was  applied  to  the 
temperature  gradients.  The  final  result  is  shown  in  Fig.  10  where  the  asymptotic  values  of  the 
shea-  band  width  2a  is  plotted  versus  logarithm  of  the  nominal  strain  rate,  log  fn.  Tne  points  depict 
vaiues  of  the  banc  w  idth  ob:ai.ned  by  deformation  gradients,  stars  are  obtained  by  the  temperature 
gradient  method.  Differences  ber.veen  those  two  approaches  increase  when  the  nominal  strain  rate 
decreases.  Tnis  is  due  to  a  m.o.'r  intensive  effects  of  heat  conduction  at  lower  strain  rates.  The  shear 
band  widths  were  found  to  have  a  minimum  2a-:.-.  “  I'^O  urn  a:  fa  =  3.2x10^  s*F  Fig.  10  shows 
also  mat  at  lower  smain  rates  the  shear  bands  are  more  diffused.  .At  high  strain  rates  inertia  and  rate 
sensitivity  may  again  increase  t.ne  widtii  of  .ASB’s.  However,  the  dimensions  of  the  widths  are  in 
general  agreement  with  e.xper.mentai  obser\  ations  for  annealed  mild  steels,  including  XC18  one. 
Som.e  steels  with  a  high  >-:eid  smess  a.nd  low  strain  hardening  may  narrow  substantially  the  ASB 
widths.  For  exampie.  it  has  been  found  in  [13]  that  the  shear  band  width  is  about  20  |Jim  for  HY- 
100  low  alloy  steel  tested  a:  f-  =  1.6x103  s‘F 

Seventeen  computer  simulations  of  deforming  layer  with  thickness  h  =  2.0  mm  and  with  adiabatic 

boundary  conditions,  for  seventeen  values  of  the  nominal  strain  rate  F.,,  80  S’^  ^  fn  ^  10^  s-^ 

have  revealed  existence  of  the  strain  rate  range  wiihm  w'hich  the  development  of  the  catastrophic 

shear  is  the  most  easiest.  Tnis  region  is  around  the  nominal  strain  rate  of  10^  s’F  Also  this  studv 

* 

has  confirmed  existence  of  the  minimum  of  both,  die  instability  strain  and  the  final  localization 
strain,  reponed  earlier  in  [34]  and  in  [32].  The  minimum  of  the  .ASB  thickness  was  also  revealed  at 
the  nominal  strain  rate  f-  =  3.2x103  s*F  Evolution  of  the  critical  values  Fg,  Floc  and  2a  at 
different  nominal  strain  rates  found  from  the  numencal  simulation  has  provided  a  strong  argument 
for  exoeriment  in  this  area. 


New  experimental  technique  of  impact  shearing 

The  torsion  tests  of  thin-walled  tubular  specimens  are  very  effective  method  to  study 
advanced  plastic  deformation,  for  both,  quasi-static  and  dynamic  loadings.  In  the  literature  three 
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types  of  torsion  devices  are  reponed  :  the  rofary  fiywh^e)  raachines,  the  high  speed  hydraulic 
testing  machines  and  the  SHTB.  All  those  systems  of  loading  have  its  limitations.  For  example, 
SHTB  technique  has  a  strong  limitations  in  maximal  strain  (the  net  angle  ot  rotation),  maximum 
nominal  strain  rate  of  the  order  F  =  2x10^  s*^  (for  a  standard  specimen  length),  as  well  the 
variation  of  strain  rate  during  the  test  are  quite  high.  At  high  strain  rates  the  finite  riseume  of  the 
incident  wave  may  cause  relatively  large  specimen  deformation  during  the  absence  of  mechanical 
equilibrium  between  incident  and  transmitter  bars. 

Another  experimental  technique,  very  useful  in  material  testing,  is  the  double  shear  test 
introduced  by  Ferguson  et  al,  (1967)  and  perfected  by  Campbell  and  Ferguson  (1970).  Since  one 
of  the  promising  specimen  geometries  to  study  dynamic  plasticity  and  ASB’s  is  the  Double-notch 
Shear  specimen,  it  has  been  decided  to  use  this  concept  in  shear  tesdng  at  medium  and  high  strain 
rates.  Originally,  Campbell  and  Ferguson  applied  the  loading  scheme  consisted  with  the  incident 
Hopkinson  bar  and  Oonsmitter  Hopkinson  tube  to  study  the  temperature  and  strain-rate  dependence 
of  the  yield  stress  of  a  mild  steel.  Interpretation  of  oscillograms  was  the  same  as  for  Split 
Hopkinson  Pressure  Bar  (SHPB  or  original  Kolsky  apparaiuj.).  Due  to  combination  of  a  small 
gage  length,  Ls «  0.84  mm,  that  is  the  shearing  part  of  tlie  DS  specimen,  and  a  standard  risedme  in 
the  incident  bar,  t,  =  20  ps,  relatively  advanced  plastic  defonnation  was  reached  during  the  risedme 
period.  Tlus  negative  feature  of  the  bar/tube  configuration  was  recently  discussed  in  [43]. 

It  was  decided,  the.a,  to  develop  a  new  experimental  technique  based  on  a  Modified  Double 
Shear  specimen.  The  new  experimental  technique  combines  several  posidve  factors  already 
experienced  with  the  other  setups.  This  new  technique  has  been  described  in  some  details  else¬ 
where,  Klepaczko,  [26, 43).  The  scheme  of  the  loading  setup  and  measuring  devices  arc  shown  in 
Fig.  11.  The  direct  impact  was  applied  to  load  the  MDS  specimen.  In  this  way  the  risetime  in 
specimen  loading  present  in  the  bar/tube  configuration  is  pracucally  eliminated.  The  flat-ended 
projectiles  of  different  lengins  made  of  maraging  steel  and  of  diameter  dp  ~  10  nun  are  launched 
from  an  air  gun  with  desired  velocities  Vq,  1  m/s  S  Vq  <  200  m/s.  Tne  impact  velocity  is  measured 
by  the  setup  consisting  of  three  sources  of  light,  fiber  optic  leads  1.  2.  3  and  three  independent 
photodiodes.  The  lime  interv^als  of  dork  signals  from  the  photodiodes  generated  during  passage  of 
a  projectile  are  recorded  by  two  time  counters.  The  setup  with  three  light  axes  makes  it  possible  to 
determine  acceleradon/deceleration  of  a  projectiie  just  before  impact,  thus,  an  e.xact  vaiue  of  Vo  can 
be  determined. 

Axial  displacement  L\  (t)  of  t.he  ccnirai  pan  of  the  MDS  specimen  is  measured  as  a  function 
of  lime  by  on  optical  cxiensometer  E,  acting  os  an  non-contact  displacement  gage.  The  optical 
extensometer  reacts  to  the  axial  movement  of  a  small  black  and  white  target  cemented  to  the  ceniml 
part  of  the  MDS  specimen. 

Axial  force  transmitted  by  (he  specimen  symmetric  supports  can  be  determined  as  a  funoion 
of  lime  from  the  transmitted  longitudinal  wave  et  (t)  measured  by  saain  gages  Ti,  DC  suf^dy  unit 
St  and  amplifier  At.  All  electric  signals  are  secorded  by  distal  oscilloscope  DO  and  next  stored  in 
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typcs  of  torsion  devices  are  reponed  ;  the  rotary  tlvA-hce)  machines,  th<;  high  speed  hydraulic 
testing  machines  and  the  SHTB.  All  those  systems  of  loading  have  its  limitadons.  For  example, 
SHTB  technique  has  a  strong  limitations  in  maximal  strain  (the  net  angle  ot  rotadon),  maximum 
nominal  strain  rate  of  the  order  f  =  2xl(P  s'*  (for  a  standard  specimen  length),  as  well  the 
variadon  of  strain  rate  during  the  test  are  quite  high.  At  high  strain  rates  the  finite  risedme  of  the 
incident  wave  may  cause  relatively  large  specimen  deformadon  during  the  absence  of  mechanical 
equilibnum  between  incident  and  transmitter  ban. 

Another  experimental  technique,  very  useful  in  material  testing,  is  the  double  shear  test 
introduced  by  Ferguson  et  al.  (1967)  and  perfected  by  Campbell  and  Ferguson  (1970).  Since  one 
of  the  promising  specimen  geometries  to  study  dynamic  plasdcity  and  ASB’s  is  the  Double-notch 
Shear  specimen,  it  has  been  decided  to  use  this  concept  in  shear  tesdng  at  medium  and  high  strain 
rates.  Originally,  Campbell  and  Ferguson  applied  the  loading  scheme  consisted  with  the  incident 
Hopkinson  bar  and  transmitter  Hopkinson  tube  to  study  the  temperature  and  strain-rate  dependence 
of  the  yield  stress  of  a  mild  steel.  Interpretation  of  oscillograms  was  the  same  as  for  Split 
Hopkinson  Pressure  Bar  (SHPB  or  original  Kolsky  apparaiui).  Due  lo  combination  of  a  small 
guge  length,  1.)  =  0.84  mm,  that  is  the  shearing  pan  of  the  DS  specimen,  and  a  standard  risedme  in 
the  incident  bar,  t,  •  20  ps.  relatively  advanced  plastic  deformation  was  reached  during  the  risedme 
period.  This  negative  feature  of  the  barAube  configuration  was  recently  discussed  in  143). 

It  was  decided,  the.i.  to  develop  a  new  experimental  technique  based  on  a  Modified  Double 
Shear  specimen.  The  new  experimental  technique  combines  several  posidve  factors  already 
experienced  with  the  other  setups.  This  new  technique  has  been  described  in  some  details  else¬ 
where,  IClepaczko,  [26. 4.3].  The  scheme  of  the  loading  setup  and  measuring  devices  axe  shown  in 
Fig.  11.  The  direct  impact  was  applied  to  load  the  .VIDS  specimen.  In  this  way  the  risedme  in 
specimen  loading  present  in  the  bar/tube  configuradon  is  pracdcally  eliminated.  The  flat-ended 
projectiles  of  different  lengins  made  of  maraging  steel  and  of  diameter  d^  >  10  rrun  are  launched 
from  an  air  gun  with  desired  velocities  V„,  1  m/s  S  S  200  m/s.  The  impact  velocity  is  measured 
by  the  setup  consisting  of  three  sources  of  light,  fiber  optic  leads  1.  2.  3  and  three  inckpendeni 
photodiodes.  The  time  intervals  of  dark  signals  fftmi  the  photodiodes  generated  during  passage  of 
a  projecdle  are  recorded  by  two  time  counters.  The  setup  with  three  light  axes  makes  it  possible  to 
determine  accekradonAlecekradon  of  a  projecdle  just  before  impact,  thus,  an  exact  value  of  V,  can 
be  deteimined. 

Axial  displacement  t\  (t)  of  t.he  central  pan  of  the  MDS  specimen  is  measured  as  a  funedoo 
of  dme  by  an  optical  extensometer  E.  acting  as  an  non-contact  displacement  gage.  The  optical 
extensometer  reacts  to  the  axial  movement  of  a  small  black  and  white  orget  cemented  to  the  cenind 
pan  of  the  MDS  specimen. 

Axial  force  transmitted  by  the  specimen  symmetric  supports  can  be  detenraned  as  a  fuoctkM 
of  time  from  the  transmitted  kmgitudinat  wave  ct  (0  measured  by  strain  gages  Ti,  DC  siqiply  unit 
S|  and  an^lifkr  A).  All  ekctric  signals  are  recorded  by  digital  oscilloscope  DO  and  next  stored  in 
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the  computer  hard  disc.  In  addinon,  a  hard  copy  of  the  recorded  signals  can  be  produced  with  an 
XY  recorder  or  a  printer.  This  new  configuration  of  experimental  seaip  permits  for  a  wide  variation 
of  the  nominal  strain  rate.  Direct  determination  of  the  axial  displacement  permits  also  for  a  more 
exact  evaluadon  of  deformadon  history,  a  very  important  piece  of  informadon  in  development  of 
adiabadc  shear  bands.  After  eliminadon  of  dme  a  force-displacement  curve  can  be  constructed  for 
the  MDS  specimen  and  next  tCT)  or  nf)  curves  can  be  determined.  The  experimental  technique 
described  above  has  appeared  quite  effeedve  and  flexible  for  material  tesdng  at  high  strain  rates  as 
well  as  for  studies  of  adiabadc  shear  bands.  Detailed  description  of  this  technique  including 
caiibradon  of  the  MDS  specimen  by  FE  is  given  elsewhere,  [43]. 

The  net  displacement  5*  imposed  on  the  central  pan  of  the  MDS  specimen  is 


5j(0  =  5^(t)-C 


c-j-  d})  di3 


(37) 


■.vhere  5^  li)  is  dispiacement  measured  by  the  optical  extensometer/displacement  gage  shown  by  E 
in  Fig.  1 1.  Ej  (t)  is  the  transmitted  wave  signal,  measured  by  the  snain  gage  T'.  cemented  on  the 
c-ansmitter  tube  sunact.  Q  is  me  elastic  wave  speed  in  the  tube.  If  h  is  the  deforming  gage  length 
of  the  MDS  specimen,  the  normnai  shear  strain  can  be  found  from  me  formula 


The  nominal  smur.  rate  can  be  found  also 


(3S) 


(59) 


It  1$  assurr.ed  in  the  elementary  approach  that  the  rariSmiaed  stgnai  17  ttl  u  dte  Hopkinson  tube  is 


sutTicient  to  rind  the  axtai  fome  tnnsmined  by  the  MDS  supports  into  die  rube 


F  i  t  ) »  £  .A 


t  E«“<t) 


(40) 


where  .Ag  is  the  cross  stetson  of  the  tube  and  £  ;s  Young’s  moduius.  The  cross  secdon  Ag  is 
related  to  the  external  and  mtemal  tube  diameter 


Ac  =  T‘D*d^ 


Tne  shear  stress  ;n  uhe  MDS  specixren  can  be  found  from  die  fcrotula, 


t(t)« 


^  Fit) 
‘*:a6 


(41) 


where  a  b  is  the  aedve  cross  secdon  of  one  suk  of  the  MDS  specimen.  ^  is  the  calibration  factor.  ^ 
2  1.  and  F  is  the  axial  force  in  the  Hopkinson  tube.  The  final  fonnula  for  the  shear  stress  is 
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T(t)  = 


71 


D‘-d^^E 


Sab 


EtO) 


(42) 


Thus,  the  shear  stress  in  the  MDS  specimen  is  proponional  to  the  current  signal  of  the 
transmitted  longitudinal  wave  ej.  All  those  relations  have  been  derived  using  the  elementary  wave 
theory  of  longitudinal  elastic  waves.  In  addition,  it  was  assumed  in  derivation  of  eq.  (15)  that  the 
local  pressures  applied  by  MDS  specimen  at  the  end  of  the  Hopkinson  tube  are  transmitted 
uniformly  at  the  distance  x  «  3D,  this  is  dynamic  St.  Venant  rule.  At  the  impact  velocity  V  =  100 
m/s  time  to  fracture  of  the  MDS  specimen  is  relatively  shon,  for  example,  for  Ter  =  1.0  the  time  to 
fracture  is  ter  ••  20  p.s,  it  gives  the  wavelength  k  ■»  100  mm  and  the  relative  wavelength  X  = 


2X/(D-d) ;  is  X  =  11.1  ;  D  =  32  mm  and  d  =  14  mm.  Those  wavelengths  indicate  that  some 

geometric  dispersion  may  be  present,  but  the  so-called  Pochhammer-Chrec  vibrations,  and  of 
course  geometric  dispersion,  are  much  smaller  in  a  tube  of  the  same  external  diameter  than  in  a 
solid  bar,  (44).  A  more  exact  analyses  of  wave  dispersion  in  tubes  are  in  progress  and  some  of 
them  are  discussed  in  the  next  pan  of  this  repwt. 

Another  imponant  parameter  is  the  length  of  projectile  Lp.  It  is  desirable  to  reach  the  time  to 
fracture  u  during  the  first  cycle  of  clastic  wave  propagation  in  a  proicctile.  If  the  nominal  suain  rate 
r  is  assumed  and  the  critical  shear  strain  Ter  can  be  estimated,  the  length  of  the  projectile  which 
assures  the  conuct  with  the  MDS  specimen  up  to  fracture  can  be  found  from  the  formula 


t  Isr 

Lp*  , 

(43) 

For  example,  when  fer »  l.O  and  F  *  5x10*  s  *  the  projectile  length  is  reduced  to  Lp  *  50 
mm.  but  for  r  =  lO^  $  *  the  length  should  be  Lp  -  2500  nom.  It  appears  that  at  lower  strain  rates  it 
is  almost  impossible  to  deform  and  fracture  the  .N£DS  specimen  during  the  contact  time  ip  » IL^C^. 
The  longest  projectile  used  in  this  study  was  Lp  »  400  mm.  For  example,  in  Fig.  12  and  Fig.  13 
are  respectively  shown  miginal  oscillogiamrocs  for  a  mild  (XC18)  and  a  hard  steel  (V AR4340). 
The  period  of  time  to  fracture  die  tnikl  steel  specimen  at  impact  \nlocit>’  V*  ■  35  mft  was  -  260  ps. 
whereas  a  more  brittle  specimen  is  broken  at  the  impact  velocity  -  160  in/^  within  •  20  pm. 

Severtheiess,  the  elementary  w-ave  theory  can  be  applied,  as  the  fost  a}^m)adi,  to  obtain  a 
complete  information  on  the  w  hole  process  of  deformation  and  fracture  of  the  MDS  specimen,  the 
following  dau  can  be  obcaiixd  from  one  test :  t  (i) ;  F  (t) ;  F  tt) ;  r  (P;  a«d  F  (F).  Since  a  two- 
channel  fast  difitai  oscilloscope  was  used  (sampling  race  IGHz).  dte  useful  potikm  of  the  digital 
memory  could  be  taken  for  analysis  ic  detetmine  all  needed  futictions.  A  computer  prognm  was 
piepaicd  »  calculate  aU  those  functions. 
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Propagation  of  elastic  longitudinal  waves  in  tubes,  the  Hopkinson  tube,  analytic 
solutions 


It  was  shown  in  the  pre%-ious  Chapter  that  for  the  case  of  a  hard  MDS  specimen  and  at  rela¬ 
tively  high  impact  velocity  the  total  time  to  fracture  may  be  as  short  as  ~  20  p-s.  When  fracture 
occurs  in  such  a  short  ame  i:  may  expected  that  some  dispersive  effects  can  be  present  during 
transmission  of  the  elastic  signals  between  MDS  specimen  and  the  strain  gage  SR,  Fig.  11.  lliis 
distance  in  the  configuration  used  was  AXg  =  170  mm,  i.e.  about  five  external  diameteis  of  the 
Hopkinson  tube. 

It  is  then  of  interest  to  snidy  the  elastic  wave  dispersion  in  tubes.  The  problem  has  been 
divided  into  tw'o  separate  studies.  The  first  one  concentrated  on  analytic  solutions  for  two  specific 
boundary  value  problems,  that  is  for  the  mixed  problem  of  imposed  pressure  on  a  semi-inrlnite 
circular  tube  and  the  next  one  is  devoted  to  the  longitudinal  impact  (initial  conditions  in  velocity). 

The  problem  of  imposed  pressure  at  the  end  of  semi-infinite  circular  tube  has  been  solved 
analytically  by  use  of  the  double  integral  transforms  (Fourier  and  Laplace  transform).  This  method 
leads  to  solutions  given  in  the  terms  of  Fourier  and  Laplace  integrals. 

A  similar  mixed-pressure  problems  of  semi-infmite  circular  rods  were  solved  previously  in 

[45j. 

The  second  problem  solved  is  the  impact  on  a  semi-infinite  circular  elastic  tube  where  the 
axial  velocity  is  imposed  instantaneously  at  the  end  of  the  tube.  The  solution  is  obtained  under 
assumption  that  the  end  of  the  tube  has  no  tangential  stresses  but  is  suddenly  set  into  morion  with  a 
constant  velocity  and  the  lateral  sun'aces  are  free  of  stresses.  The  analogous  problem  for  a  semi- 
infinite  bar  was  solved  in  [46],  Tnc  method  applied  in  [46]  was  to  split  the  problem  into  two  pans. 
In  ihe  first  part  the  collision  of  two  half-spaces  was  considered.  Ln  the  second  pan  an  infinite  bar 
was  considered  with  an  equal  and  opposite  stresses  travelling  along  the  lateral  surface.  Using  the 
principle  of  superposition  a  stress  free  lateral  surface  was  ensured.  The  method  presented  in  [46] 
has  been  adopted  to  the  problem  of  a  semi-infinite  circular  mbe  loaded  by  impact. 

Both  solutions  for  a  tube,  which  are  mentioned  above,  are  reponed  in  details  in  the  Appendix 
N°  2  to  this  Repon,  [47]. 

Here  only  a  brief  description  of  those  solutions  will  be  given. 

A  semi-infinite  cvlindrical  elastic  rube  of  inner  radius  a  and  outer  radius  b,  which  are  referred 

„  v// 

to  the  cylindrical  coordinates  (r,  9,  z),  where  the  z-axis  coincides  with  thc/)f  the  tube,  is  conside¬ 
red.  Let  X,  )4  be  Lame  constants  and  p  be  the  mass  density  of  the  tube  material.  The  motion  is 
restricted  to  the  axisymmetric  case,  thus  the  displacement  vector  12  is  independent  of  the  circum¬ 
ferential  variable  0  and 


12  (r,z,t)  =  (Ur  (r.z,t)) ;  0  ;  U2  (r,z.t) 


(44) 
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The  equations  of  motion  of  the  linear  elasticity  in  the  cylindrical  coordinates  take  the  form 


(x  +  2n) 


I  2  2 

9  u.  1  9ur  u,  9  u 
+ - - — -  + 


r  3r  9r9z 


9  Uj  9  Uj 


\Bz 


2  9r9z 


d\ 


I . 


(X  +  2n) 


^9r9z  ^  9z  Bz^j 


1 

-  II  — . 

9uj 

-ii) 

9z 

\ 

-srj 

U2  .2  ^ 

d  Uf  d  Uj 


9z  9r  ar* 


=  P 


9t‘ 


(45) 


Consider  the  boundary  conditions 


Oil  (r.z=0.t)  =  -  f  (t) . 


and 


Ur  (r.z=0.t)  =  0 ,  for  r  e  <a.b>  .  t  e  (O.oo) 


Of.  tr=a.2.tl  =  Orr  tr=b.z.t)  *  0  . 


o-j  (r=a.2.:)  =  <Jr,  (r=b.z.t)  =  0 .  for  z  e  tO.—) .  t  €  tO.*) . 


(46) 


(47) 


where  f  t.)  is  a  given  tuncnon.  provided  that  f  ii)  «  0  for  i  <  0.  A  -urther  condition  is  that  the 
displacements  and  the  soesses  vanish  at  infinity,  a  -*  «».  Assuming  that  the  semi-infinite  circular 
tube  is  at  rest  prior  to  nme  instant  t 0.  equations  (46)  and  (47)  are  supplenxnted  by  the  initial 
conditions 


oq 

u.  ir.z.t«0)  ■-rr  tr.i-t»0)  ■  0 . 

<7\ 

(AS) 

u  j  (r.z,t«0)  ■  -rr  «  0  .  for  r  €  <a.b>  .  z  «  (0,»*) . 


The  set  of  equations  (45)  along  with  die  initial  and  boundary  cooditiotu  (46)  and  (47)  has  been 
solved  by  introduction  of  die  displaceaaent  potentials,  details  are  given  in  [46].  The  solution  by 
integral  nnsforms  has  been  found  using  the  inversion  theorrms  of  the  Laplace  andFootier  aaAs- 
fixins.  the  solution  can  be  wTiaeo  in  the  form 


Gzj(r,^p)  e^su:  (^)  dp  dc . 


The  function  f  (.),  if  the  faouncsrs'  condition  is  given  as  (46),  is  yet  undefined.  The  possibility  of 
the  followuig  cases  of  impulses  can  be  mentioned  ; 

a.  infinitely  shon  impulse,  f  r.)  =  •  fo  5  (t),  where  f©,  f©  >  0,  is  a  constant  and  5  (.)  is  the  Dirac 
function  : 

b.  infinitely  long  impulse,  f  (ty  =  -  f©  H  (t),  where  H  (.)  is  the  Heaviside  step  function  ; 

c.  impulse  of  finite  length,  f  =  -  f©  H  (t  - 1),  where  t,  T  >  0,  is  the  duration  of  the  impulse. 

The  second  boundar>-  value  problem  studied  was  the  longitudinal  impact  which  was  reduced 
to  the  mixed  condidons  of  shock  pressure  at  the  end  of  the  semi-infinite  circular  elastic  tube.  Two 
semi-infinite  circular  tubes,  moving  in  opposite  directions  with  velocity  v,  are  assumed  to  contact 
at  the  instant  t  =  0  and  at  the  plane  2  =  0.  The  lateral  surfaces  of  tubes  are  free  of  stresses. 

The  solution  of  the  problem  considered  is  devided  by  appUcadon  of  the  superposition  princi¬ 
ple  for  two  separate  problems.  The  first  of  them  is  limited  to  two  elastic  semi-spaces,  moving  in 
opposite  directions  with  velocity  v.  They  are  assumed  to  contact  at  the  instant  t  =  0  and  at  the  plane 
2  =  0.  After  time  t  =  0  the  semi-spaces  are  assumed  to  behave  as  a  single,  solid  space  which  are 
assumed  to  behave  for  t  >  0  as  a  single,  infinite  tube.  The  solution  of  this  problem  is  given  in  [46] . 
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Ur  (r,2,t)  =  0  , 


Uz  (r,z.t)  = 


-i-  V  t  for  ±  z  >  c  ti 

I-  V  z/C)  for  I  z  1  <  cit 


The  displacements  assumed  above  induce  the  following  stresses,  [47] 
azr  (r,z,t)  =  0  , 

fO  for  I  z  I  >  cit 


cJzz  (r,z,:)  = 


-  (}i-i-2|.t)  v/ci  for  I  z  I  <  cit 


(r,z,t)  = 


rO 


-  A  V/Cl 


for  !  z  I  >  cit 
for  !  z  I  <  Cit 


(50) 


(51) 


It  is  seen  trom  equations  (51)  that  the  solution  of  the  first  problem  causes  non-zero  radial  tractions 
on  the  lateral  surfaces  of  the  circular  tube.  Applying  the  principle  of  superposition  to  ensure  a 
stress  free  lateral  surfaces  of  the  infinite  tube  the  second  problem  is  considered,  which  is  defined 
by  equations  of  morion  (45).  the  boundary  conditions 


azr  (r=a.z,t)  =  Ozr  (r=b,z,t)  =  0  , 

I® 

Orr  (r=a.z,t)  =  On  (r=b,z,t)  =  < 

i-  X  v/ci 

for  z  €  (-<»,-:-<»),  t  €  (0,  +  ®«) , 


for  I  z  I  >  Cit 
for  I  z  I  <  cit 


(52) 


and  the  initial  conditions  given  in  (5c)  for  z  e  (- »«,  +  ®«).  A  further  condition  is  that  the  displace¬ 
ments  and  the  stresses  vanish  when  I  z  ! 

The  solution  of  the  problem  is  found  by  using  displacement  potendaiy  and  Fourier  exponen¬ 
tial  iransfonns,  as  well  as  Laplace  transform.  Using  the  inversion  theorems  of  the  Laplace  and 
Fourier  exponential  transforms  the  solution  can  be  written  in  the  form 


u,(r.z.t)  =  - 


4?:  i 


_F 


'Z  Dt 


u.  (r,c.p)e  •  e'  dp  , 


u,(r,z,t)  =  - 


47:  i 


Uj  (r,^p)e''^e“dp  dq. 


J  . 


a^(r,2,t)  =  - 


a^(r4p)e'^e^dp  d^, 


4::  i 


(53) 


G^(r,z.t)  =  — 
4;c  i 


Gjr(r,^p)e‘'^e^dpdq, 


G_(r,z.:)  = 


_  .  -  .  cz  ot  , 
a_(r.c.p)e  '  c  dp  dc , 


The  soluiions  (44)  and  (53)  enable  to  calculate  effects  of  dispersion  for  idealised  impulses,  type 
Dirac  and  square  of  infinite  and  finite  lengths. 


Propagation  of  elastic  longitudinal  waves  in  tubes,  the  Hopkinson  tube,  numerical 
study 

Although  analytic  solutions  have  been  found  for  an  elastic  semi-infinite  tube  loaded  instanta¬ 
neously  by  pressure  or  by  velocity,  the  real  case  with  MDS  specimen  attached  to  the  front  of  the 
Hopldnson  tube  is  still  more  complicated,  mostly  due  to  non-symmetric  circumferential  distribution 
of  pressures.  Such  a  case  is  impossible  to  solve  analytically,  and  hence,  scries  of  numerical  studies 
have  been  undertaken  using  a  dynamic  FE  code.  A  non-symmetric  distribution  of  pressures  have 
been  assumed  at  the  face  of  the  tube  in  the  form  of  two  contact  areas  of  the  MDS  specimen,  that  is 
a  =  R  ■  r,  a  =  9.0  mm,  where  R  is  the  outer  radius  and  r  is  the  inner  radius,  b  =  6.0  mm,  where  b 
is  the  thickness  of  the  MDS  specimen.  The  total  area  by  which  the  total  force  is  transmitted  is  A  « 
2ab.  This  is  shown  in  Fig.  14.  In  Fig.  15  are  shown  all  cases  of  trapezoidal  pulses  assumed  at  the 
end  of  the  tube  on  the  shadowed  surfaces  2ab. 
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The  following  material  constants  were  assumed  for  the  Hopkinson  tube,  E  =  2.1  *  10^  MPa, 
V  =  0.3  and  p  =  7.85  g/cm^.  The  Hopkinson  tube  has  been  modelled  as  semi-infinite  by  its 
anificial  extension  with  a  very  low  wave  speed  and  the  same  impedance.  Thus,  all  reflections  from 
the  other  side  of  the  tube  have  been  eliminated.  Comparison  of  calculations  with  a  normally  exten¬ 
ded  tube  by  FE  (the  same  Co,  and  p)  revealed  only  small  differences  between  those  two  cases. 

Since  the  main  interest  is  concentrated  on  high  impact  velocities  imposed  on  MDS  speci¬ 
mens,  Vq  =  100  m/s,  the  cases  N°  1  and  N°  4  will  be  briefly  discussed  in  this  Report.  Case  N®  1  is 
assumed  as  pmax  =  450  MPa,  tr  =  20  |is  ;  to  =  10  ps  and  t,  =  20  ps.  The  results  of  FE  calculations 
in  the  form  of  longitudinal  strain  ei(t)x  ande2(t)x,  where  fii  is  the  longitudinal  strain  at  the  outer 
surface  of  the  Hopkinson  tube,  and  £2  at  the  inner  surface,  are  shown  in  Fig.  16a/16b  for  xi  = 
90.5  mm  (-  3D  of  the  tube)  and  in  Fig.  17a/17b  for  X2  =  170.5  mm  (~  5D).  In  each  figure  a  and  b 
the  results  are  shown  for  the  external  surface  of  the  tube  (solid  line)  and  the  internal  surface 
(broken  line).  In  addition,  figures  a  are  obtained  for  the  vertical  plane  from  Fig.  14a,  that  is  the 
plane  in  which  the  MDS  specimen  is  attached.  Whereas  figures  b  refer  to  the  time-distribution  of 
longitudinal  strains  £i  and  £2  at  the  same  distances  xi  and  X2  but  with  90°  orientation. 

Those  figures  indicate  the  existence  of  higher  modes  of  vibrations  excited  by  a  non-symme- 
tric  impact.  This  iS  visible  for  both  distances  studied.  The  FE  calculations  reveal  also  small  diffe¬ 
rences  between  dme-distributions  of  longitudinal  strain  obtained  on  external  and  internal  surface  of 
the  Hopkinson  tube.  The  origin  of  those  differences  are  unknown,  and  they  may  be  caused  by  the 
FE  technique  itself.  Funher  studies  are  needed  at  this  point,  perhaps  with  different  meshing.  The 
assumed  risetime  uf  the  incident  pulse  was  tr  =  20  ps.  this  risetime  is  slightly  disioned  at  the 
beginning  for  both  distances  xj  =  90.5  ram  and  xi  =  170  mra  by  a  small  inflection.  The  descen¬ 
ding  pan  of  the  incident  pulse,  assumed  originally  as  c,  s  20  ps,  is  slighdy  delayed.  The  results 
obtained  for  the  distance  x;  s  17C.5  mm  do  not  differ  much  from  the  previous  ones  for  xi.  The 
transmitted  wave  determined  from  the  tube  surface  is  equally  dispersed,  and  the  differences  in 
waves  between  exieraai  and  internal  surfaces  are  small.  Again,  this  needs  funher  studies  for  less 
idealized  impulses. 

The  most  imponant  conclusion  drawn  from  the  FE  calculation  of  the  case  N°  1  is  that  the 
level  of  V.  e  mean  stress  of  the  plateau  determined  from  the  surface  of  the  Hopkinstm  mbe  agrees 
with  the  assumed  level  of  pressure  pm«x  multiplied  by  the  ratio  of  suna'^s. 

TTb  (54) 

The  ratio  4  in  the  real  case  of  Hq)kinson  mbe  is  4  *  0~2.  This  fKt  peniiits  for  detenmnation  of  the 
critical  stress  imposed  oa  the  .MDS  specimen  via  the  peak  of  strain  measured  by  the  Hopkioson 
mbe.  in  the  real  case  the  strain  gages  are  cemented  at  x  •  170  mm. 

The  most  interesting  case  is  the  N°  4,  because  (luring  testing  of  hard  steels  the  order  of  tr,  to 
and  4  are  similar  as  assumed  for  the  N”  4 : 4  »  5  ps.  t^  «  5  ps  and  >  5  ps.  however,  pnc\  ~  4S0 
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MPa  is  rather  a  small  value  in  comparison  to  the  real  ones  for  a  hard  steels  with  the  flow  stress  a  = 
1000  MPa.  The  results  of  those  FE  calculations  arc  shown  in  Fig.  18  and  Fig.  19.  The  same 
approach  was  used  as  tor  the  case  N°  1.  .Again,  the  auialysis  of  amplitudes  indicates  that  even  for 
such  shon  pulses  the  ma.'dmum  stresses  determined  in  MDS  specimen  after  the  back  extrapolation 
are  correct.  At  the  distance  xi  =  90.5  mm  the  rising  part  of  the  pulse  shows  at  the  beginning  a 
small  inflection,  overall  this  figure  resembles  the  oscillogram  for  VAR4340  steel  shown  in  Fig.  13. 
■At  the  distance  x:  =  170.5  mm  the  pulse  is  slightly  wider  due  to  geometry  dispersion,  but  the 
mean  value  of  strain  also  permits  for  back  extrapolation  to  find  MDS  stresses.  The  tails  of  the 
pulses  for  both  distances  show  large  amplitude  vibrations  which  resemble  to  some  extend 
measurements  (Fig.  13).  Since  the  assumed  values  of  pmax  were  relatively  low  in  comparison  to 
the  real  levels  of  stresses  in  tested  4340  steel,  another  series  FE  calculations  have  been  run  for 
Cases  N’  1  and  4.  but  with  =  10(X)  MPa  and  the  same  values  of  t.,  to  and  tj.  Those 
calculations  have  confLnned  validity  of  the  procedure  of  back  extrapolation  of  stresses  for  MDS 
specimens.  The  relative  level  of  vibrations  is  reduced  in  comparison  to  the  main  pulse.  This  is  also 
obser.'ed  for  the  tails  of  the  calculated  waves.  Other  calculations  gave  vcr>'  similar  results. 

The  main  conclusion  drawn  from  FE  dynamic  calculations  is  the  confirmation  of  the  back 
extrapolation  procedure  to  find  saesses  imposed  on  MDS  specimens  by  direct  impact  using  eq. 
^54). 


.Numerical  analysis  of  Modified  Double-Shear  specimen  and  specimen  calibration 

Tne  new  configuranon  of  expetiment  with  the  A-IDS  specimen,  as  it  is  shown  in  Fig.  12  has 
been  developed  in  LP.MM-Metz.  The  modified  specimen  geometries  are  shown  in  Fig.  20a  and 
20b.  The  modified  geometries  have  a  wider  gage  length,  U  =  2.0  mm.  in  comparison  to  [48], 
which  substantially  reduce  stress  concentrations  near  the  notch  comers. 

The  external  part  are  longer  and  those  pans  are  clamped  to  the  transmitter  tube  to  prevent 
rotations.  The  MDS  geometry  was  analysed  by  FE  technique  at  different  stages  of  plastic  deforma¬ 
tion.  The  total  number  of  elements  for  half  of  the  MDS  specimen  was  640.  Around  and  within 
deforaiing  pan  of  the  specimen  the  mesh  was  much  denser,  460  elements.  The  central  pan  of  the 
MDS  specimen  was  assument  to  depiacc  uniformly.  Since  the  width  of  the  plastically  deformed 
pans  of  the  ATDS  specimen  is  only  6.0  mm  and  the  shear  stress  z,  is  supposed  to  be  determined  for 
an  infinite  layer,  a  calibration  factor  q  must  be  introduced,  thus 

^  =  (55) 

where  z^  is  the  mean  stress  in  the  deformed  section  of  the  MDS  specimen.  Distribution  of  the  shear 
stresses  at  the  nominal  shear  strain  Ts  =  2.5  x  lO'^  is  shown  in  Fig.  21a  whereas  Fig.  21b  shows 
distribution  for  Fs  =  0.5.  Distributions  of  the  shear  strains  for  the  same  values  of  Fj  are  shown  in 
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Fig.  22a  and  22b.  All  FE  calculations  have  been  performed  for  the  quasi-static  stress-strain 
relation,  which  depicts  more  or  less  exactly  tensile  behavior  of  low  alloy  steel  XC18  at  strain  rate  e 
=  10'^  s‘i 


a  =  Gy  +  B 


I 

E 


(56) 


with  the  yield  stress  and  strain  :  Oy  =  240  MPa,  £y  =  1. 165  x  10’^ ; 


D  = 


1 


n  1  -n 
E 


P  = 


1 

1  -  n 


n  =  0.3  :  E  =  2.06  x  ICP  MPa  ;  B  =  451.1  MPa.  This  is  very  useful  form  of  a(e)  relation  since  it 
has  a  smooth  transition  from  the  elastic  to  plastic  regions,  that  is  the  first  derivative  0  =  do/de  has 
no  jumps.  Of  course,  in  the  elastic  region  a  =  Ee  for  e  <  Ey.  The  mean  value  of  the  stress  calibra¬ 
tion  factor  ^  obtained  after  FE  analysis  for  MDS  specimen  shown  in  Fig.  20a  is  ^  1.2.  It  was 

found  that  ^  slightly  increases  from  q  *»  1.12  at  F,  =  0.05  to  q  •»  1.2  at  F,  =  0.50,  so  the  mean 
value  for  0.05  S  Fj  ^  1.0  is  ~  1.2.  The  distribudons  of  shear  strain  shown  in  Fig.  22a  and  22b 
demonstrate  that  the  strain  field  is  quite  uniform  with  a  small  asymmeoy  due  to  bending  component 
from  the  supports.  It  was  also  found,  in  the  real  tests  as  well  as  from  the  FE  analysis,  that  the 
exterior  pans  of  the  specimen  which  are  supponed  also  deform  plasdcally,  causing  an  increasing 
tendency  for  rotation  (MDS  specimen.  Fig.  20a).  This  is  why  the  specimen  geometry  has  been 
slightly  modified  as  it  is  shown  in  Fig.  20b.  The  width  of  the  suppons  is  enlarged  to  10  nun  and 
the  notch  is  slightly  diper.  It  is  clear  that  due  to  modifications  introduced  into  the  specimen  geome¬ 
try  the  axiai  force  versus  a.xial  displacement  provides  a  reliable  input  data  to  find  quasi-stadc  or 
dynamic  x(F)  response. 


Preliminary  experiments,  mild  steel 

It  has  been  decided  to  perform  preliminary*  series  of  tests  on  die  mild  steel  in  the  recrystaUi- 
sed  state.  This  steel  has  been  assumed  as  a  model  material  for  the  numerical  calculations  and  all 
mechanical  tests.  In  order  to  perform  quasi-stadc  tests  and  also  tests  with  higher  strain  rates,  of  the 
order  10^  s  ^  a  special  loading  rig  has  b^n  constructed.  The  rig,  shown  in  Fig.  23.  having  its 
own  measuring  systems  (force  is  measured  by  an  SR-gages  and  signal  condidooer.  and  displace¬ 
ment  by  two  LVDT-s  and  amplifiers),  can  be  fitted  to  any  standard  tesdng  machine.  The  signals 
F(t)  and  6i(t),  could  be  stored  in  a  digital  memory,  where  F  is  the  axial  force,  8i  and  $2  ue 
LVDT  displacements.  Shear  soess  and  shear  strain  can  be  obtained  from  the  fonnulas 
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F(t) 


2  a  b 


and 


(57) 


where  a  and  b  are  respectively  the  width  of  the  defora:ed  pan  and  specimen  thickness,  the  mean 
axial  displacement  is 


5s(t)  = 


5i  (0  +  62(0 


and  r(0  = 


5s  (0 


(58) 


where  6^  is  the  mean  axial  displacement.  Elimination  of  time  from  t(t)  and  r(t)  pennits  for  deter- 
minadon  of  ^(nf-  In  addition,  the  strain  rate  must  be  determined  from  the  formula 


r(t)  = 


1  d5;(t) 
Is  d  t 


(59) 


.Application  of  the  closed-loop  dynamic  tesdng  machine  makes  it  possible  to  reach  for  the  MDS 
specimen  the  strain  rate  of  the  order  F  •  10^  S'^.  Of  course,  the  maximum  strain  rate  is  a  linear 
function  of  the  axial  velocity  and  proportional  to  I/l*.  Fig.  24  shows  log  f  as  a  function  of  gage 
length  1,  for  two  axial  velocides :  20  m/s  and  100  m/'s.  By  the  broken  venical  lines  arc  marked  the 
gage  lengths  1$  =  0.25  mm  and  0.84  tim,  used  in  (48,  49],  and  I,  =  2.0  mm  used  in  the  MDS  spe- 
ctmen.  For  example,  in  the  case  0:  impact  velocity  100  m  s  and  MDS  specimen  the  ncminal  strain 
rate  F  is  5  x  l(y  s*‘  is  reached.  Thus,  the  new  experimental  technique  can  cover  also  the  high  strain 
rate  region.  10^  s'*  ^  F  S  10^  s’*  (2.0  m/s  ^  V  <  200  m/s). 

First  scries  of  experiments  were  performed  on  low-alloy  steel  .XC18  (Trench  Standards)  of 
the  following  compos’^on  in  ^  :  C  0.17  ;  Mn  0.5S  ;  Si  0.21  ;  S  0.032  ;  P  0.024.  The  specimen 
material  was  supplied  in  the  form  of  cold  drawn  cars  DI.A  12  mm.  .After  precision  machining,  all 
specuneas  were  vacuum  an.-tealed  a;  -  1000  K  tc:  I*  -  hrs  and  furnace  cooled.  Meiallographic  exa- 
iTunation  showed  a  uniform  microscructure  w;th  evenly  distributed  grains  of  ferrite.  The  tests  on 
.MDS  specimens  were  performed  withui  two  ranges  of  strain  rates :  medium  rates  10^  S’^  S  F  S10^‘ 
s  *  with  closed-loop  fast  tcstmg  machine  and  high  rates  Ky  S’‘  S  F  S  iO^  s’‘  with  the  new 
experimental  technique.  One  of  oscillograms  obtained  with  impact  velocity  V  •  •  37  tn/s  is  shown 
i.*;  Fig.  12a.  The  upper  trace  cotrss  from  the  srain  gages  on  the  ransmit^r  tube  and  the  lower  one 
from  me  optical  extensoineter.  It  con  be  seen  from  F:g.  12a  that  afarr  the  upper  yield  peak  there  is  a 
small  oscillation  on  the  upper  trace,  they  are  believed  to  be  panly  caused  by  the  out-of-phase  longi¬ 
tudinal  vibrations  of  the  transmitter  tube  revealed  by  the  FE  dyrusaic  calculadons.  Since  the  length 
of  the  projectile  was  «  180  mm  with  the  droe  of  contact  tp  -  74  ^s.  armind  three  fuU  cycles  of 
wave  reflexions  had  occured  during  specimen  deformation.  They  are  recognizable  in  the 
oscillogram.  As  a  whole,  both  the  upper  and  lower  yield  stresses  can  be  found,  as  well  u  the 
piasnc  flow  curve  up  the  fracture  by  adiabanc  sheanng.  The  signal  from  the  optical  extensometer 
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showfthat  the  displacement  rate  of  the  central  part  of  the  MDS  specimen  is  almost  constant.  The 
optical  contact  was  lost  txsfore  the  specimen  fractured. 

A  number  of  MDS  specimens  was  tested  at  different  impact  velocities  from  30  m/s  to  about 
100  m/s.  The  upper  yield  stress  in  shear  Xy  is  plotted  in  Fig.  25  against  the  logarithm  of  shear 
strain  rate  T.  In  the  region  of  low  strain  rates,  values  of  shear  stress  Xy  were  obtained  from  tensile 
tests,  [39],  and  transformed  into  shear  by  Xy  =  Oy/-^  and  F  =  V3  £.  In  this  region  the  rate  sensiti¬ 
vity  P  is  characterized  by  a  small  nearly  constant  slope ;  P  =  dty/3log  F,  [31, 48, 50).  In  the  region 
of  high  strain  rates  a  substantial  increase  of  the  upper  yield  stress  is  observed.  This  is  shown  in  the 
linear  scale  of  F  in  Fig.  26.  This  feature  wiU  be  discussed  in  the  further  pan  of  this  Repon.  In  Fig. 
27  the  mean  values  of  the  lower  yield  stress  tyi  (vibrations  are  eliminated),  obtained  for  XC18  steel 
are  compared  with  results  for  a  mild  steel  reponed  by  Campbell  and  Ferguson  [48] ;  composition 
of  En3B  (British  Standards)  steel  as  given  in  [48]  is :  C  0.12 ;  Mn  0.62 ;  Si  0.10 ;  S  0.029.  Becau¬ 
se  of  the  higher  carbon  content  in  XC18  steel  the  level  of  the  lower  yield  stresses  is  overall  higher 
than  for  En3B  steel.  On  the  other  hand  DS  specimen  used  in  [48]  with  slots  0.84  mm  and  length  of 
deformed  pan,  a  =  3.2  mm,  should  have  a  much  higher  stress  calibration  factor  ^  in  comparison  to 
the  MDS  specimen  where  q  -  1.2.  It  is  also  probable  that  the  instantaneous  strain  rates  in  DS 
specimen  tested  in  the  SHB  (bar  and  tube)  were  much  lower  than  the  nominal  values  used  in  plot¬ 
ting  of  Fig.  27.  This  couid  be  a  reason  why  tests  with  the  new  technique  give  higher  values  of  tyi. 

Another  comparison  was  made  between  results  obtained  for  a  high  punity  iron  (99.99%  Fe) 
in  [51],  where  pressure-shear  impact  e.xperiments  were  used.  Rat  specimens  with  thickness 
between  300  and  3  pm  were  sandwiched  between  two  hard,  parallel  plates  that  are  inclined 
relative  to  their  direction  of  approach.  Nominal  stresses  and  strains  in  the  specimens  were  determi¬ 
ned  from  elastic  wave  prorlles  monitord  at  the  rear  surface  of  one  of  the  hard  plates.  It  is  seen  in 
Fig.  28.  where  the  compaiison  is  shown,  that  at  high  shearing  rates  the  shear  stresses  are  of  the 
same  order,  however,  similar  values  for  tyi  obtained  from  .MDS  specimen  occur  at  strain  rate  f  • 
10^  $*’  and  those  obtained  from  pressure  shear  e.xperiment  at  f  •  10^  s  ^  General  trends  in 
changes  of  the  lower  yield  stress  are  similar  for  both,  so  different,  experimental  methods. 

It  is  w*eU  known  that  in  the  region  of  low  and  medium  strain  rates.  10**  s  *  F  2  10^  s  ^  the 
rate  and  temperature  sensitivities  of  flow  stress  are  controlled  by  the  thermally  acuvated  kinetics  of 
dislocations.  [38, 48.  50).  In  this  region  the  rate  sensitivity  ^  « (dt/diog  Or  does  not  depend  on 
strain  rate  f .  Howev'cr.  a*  it  is  seen  from  Fig.  25  and  Fig.  27.  the  high  strain  rate  region,  typically 
10^  s*’  ^  F  £  10^  s  ^  is  characterized  by  a  rapid  increase  of  ^  with  increasing  strain  rate.  It  is 
commonly  assumed,  and  it  was  verified  in  many  theoretical  and  experin»ntal  studies,  that  in  BCC 
structures  the  instantaneous  rate  sensitivity  (3  is  the  result  of  dislocation  movement  by  the  mecha¬ 
nism  of  the  thermally  activated  formation  of  kink  pain  over  Peterls  pocntial.  for  example,  (38. 40, 
41.  52).  Since  this  approach  is  well  known,  derailed  discussion  of  its  applications  will  be  here 
omiaed. 
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When,  within  the  region  of  high  strain  rates,  the  yield  stress  or  flow  stress  are  plotted  in  a 
linear  scale  of  strain  rate,  the  experimental  data  can  be  represented  by  a  stright  line,  148].  Experi¬ 
mental  points  for  XC18  steel  from  Fig.  25  have  been  transformed  in  the  linear  scale  of  f  and  the 
final  result  is  shown  in  Fig.  29  for  the  case  of  the  upper  yield  stress  and  in  Fig.  30  for  the  lower 
yield  stress.  Indeed,  in  both  cases  the  linear  approximation  fits  experimental  points.  The  following 
linear  relation  was  analysed  by  the  least  squares  method 

(60) 

where  ty  is  the  threshold  stress  and  ri  is  the  pseudo- viscosity.  Both  yield  stresses  were  analysed 
and  the  best  linear  fit  is  shown  in  Figs.  29  and  30  by  the  broken  lines.The  complete  results  of 
those  analyses  are  given  in  Table  2,  whereas  the  reanalysed  results  from  ref.  [48]  are  given  in 
Table  3. 


Table  2 


Steel  XC18 

rt  [kPa.s] 

Upper  y.  stress 

14.3 

0.954 

Lower  V.  stress 

4.0 

0.870 

Table  3 


Steel  En  3  B 

tv  [MPa] 

n  [kPa.s] 

tt 

Lower  y.  stress 

241.4 

2.1 

0.960 

r  =  0.075 

241.7 

2.5 

0.911 

r  =  o.:i 

246.2 

2.5 

0.924 

r  =  0.38 

249.2 

2.6 

0.977 

Comparison  of  Tables  2  and  3  shows  that  values  of  for  the  lower  yield  stress  are  of  the  same 
order,  although  MDS  specimen  loaded  by  direct  Lmpact  show  r) -values  alnx)St  twice  higher.  Becau¬ 
se  of  die  direct  impact  loading,  with  the  rise  ume  I  pis  5  tj  S  2  ^s.  behaNnor  of  the  upper  yield 
smess  is  different  for  MDS  specimens.  It  is  associated  with  the  delay  time  of  plasticity*  in  mild 
steels.  (50). 

An  early  inieipreiadon  of  the  linear  relation  t  (f)r,T  is  that  the  stress  intercqits  at  f  «  To :  that 
is  i  (Or.  can  be  detemuned  by  the  temperature -dependent  threshold  stress.  When  the  applied 
stress  T  exceeds  the  tlireshold  stress.  1 2  t.  the  dislocation  velocity  is  controlled  by  dissipation  of 
energy  when  dislocation  moves  through  the  lutice.  It  is  commonly  assumed  that  this  dissipation  is 
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of  linear  viscous  nature,  and  then  the  stress  component  would  be  proportional  to  the  mean  dis¬ 
location  velocity  v  .  Then  the  viscous  drag  relation  may  be  expressed  in  the  form 


X 


n 


(61) 


where  B  is  a  damping  constant.  When  the  Orowan  relation  is  introduced  into  eq.  (61),  the  compo¬ 
nent  x,^  can  be  expressed  in  terms  of  the  mean  mobile  dislocation  density  pn,.  and  F  =:  Pnj  b  v  , 


or 


(62) 


where  t]  =  B/b^  pm  is  the  coefficient  of  viscosity. 

This  interpretation  has  many  weak  points.  The  first  one  is  that  the  thermally  activated  proces¬ 
ses  are  still  operative  at  high  strain  rates.  A  more  exact  analysis  of  experimental  data  obtained  by 
Campbel  and  Ferguson,  [48],  were  published  in  (53).  It  was  concluded  in  this  latter  analysis  that 
the  whole  region  of  strain  rates  covered  in  [48]  can  be  approximated  by  an  unique  constitutive  rela¬ 
tion  based  on  the  thermally  activated  micromechanism  of  double  kink  kinetics.  This  analysis  indi¬ 
cates  on  the  possibility  that  the  theory  of  thermally  activated  flow  is  exact  enough  to  describe  rate 
sensitivity  of  iron  and  mild  steels  up  to  strain  rates  of  the  order  10^  to  10®  s*K  [53]. 

.Another  possibility,  and  more  recent,  to  explain  the  pseudo-viscosity  lies  in  the  fact  that 
some  experimental  data  clearly  indicate  that  the  mechanical  threshold  %  is  itself  rate  sensitive,  [54, 
55],  Theoretical  analyses  of  structural  evolution  which  can  lead  to  a  rate  sensitive  mechanical 
titreshold  were  also  published  recently,  [56, 57J. 

An  imponant  observation,  directly  related  to  ASB  in  XCIS  steel,  are  values  of  the  critical 
total  strain  (elastic  plus  plastic)  detemtined  as  a  function  of  imposed  velocity  on  N£DS  specimens. 
The  final  result  is  shown  in  Fig.  31  where  the  critical  shear  strains.  Fm  -  6r,/1s,  ore  shown  as  a 
function  of  log  F.  The  mean  values  ore  shown  as  circles.  Tltc  results  for  two  lowest  strain  rates,  F 
«  75  s*^  and  f  »  750  s*^  were  obtained  with  the  loading  rig  and  the  fast  closed-loop  machine,  the 
data  for  three  highest  strain  rates  in  the  range  l.S  «  10^  S‘^  ^  f  4.7  •  10^  s*^  were  obtained  with 
the  new  experimental  technique.  The  impact  velocity  at  the  highest  nominal  strain  rate  was  - 100 
m/s. 

Since  the  dimensio.ns  of  deforming  parts  of  the  MDS  specimen  are  the  same  as  the  thickness 
of  the  layer  studied  by  numerical  medtods,  and  the  same  steel  was  assumed  for  the  numerical 
study,  a  direct  comparison  is  possible  between  experimental  results  shown  in  Fig.  31  and  numeri¬ 
cal  results  of  Fig.  3.  Although  values  of  localization  strains  are  overall  smaller  than  the  critical 
strains  found  by  experiment,  at  tower  strain  rates  the  same  trend  is  observed,  that  is  the  total  shear 
strain  increases  as  a  function  of  log  F.  However,  it  is  clear  that  at  the  highest  strain  rates  an  cdmipt 
drop  in  the  critical  strain  occurs.  For  example,  the  mean  value  of  the  critical  strain  fiom  six  tests 
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of  linear  viscous  nature,  and  then  the  stress  component  Tt,  would  be  proportional  to  the  mean  dis¬ 
location  velocity  v  .  Then  the  viscous  drag  relation  may  be  expressed  in  the  form 


where  B  is  a  damping  constant.  When  the  Orowan  relation  is  introduced  into  eq.  (61),  the  compo¬ 
nent  can  be  expressed  in  terms  of  the  mean  mobile  dislocation  density  and  F  =  P;n  b  v  , 


or  t^  =  Tir 


where  "n  =  B/b^  is  the  coefficient  of  viscosity. 

This  interpretation  has  many  weak  points.  The  first  one  is  that  the  thermally  activated  proces¬ 
ses  are  still  operadve  a:  high  strain  rates.  A  more  exact  analysis  of  experimental  data  obtained  by 
Campbel  and  Ferguson,  [48],  were  published  in  [53].  It  was  concluded  in  this  latter  analysis  that 
the  whole  region  of  strain  rates  covered  in  [48]  can  be  approximated  by  an  unique  consonidve  rela¬ 
tion  based  on  the  thermally  activated  micromechanism  of  double  kink  kinetics.  This  analysis  indi¬ 
cates  on  the  possibility  that  the  theory  of  thermally  activated  flow  is  exact  enough  to  describe  rate 
scnsidvir/  of  iron  and  mild  steels  up  to  strain  rates  of  the  order  lO®  to  10*  S'^  [53]. 

.Another  possibility,  and  more  recent,  to  explain  the  pseudo-viscosity  lies  in  the  fact  that 
some  experimental  data  clearly  indicate  that  the  mechanical  threshold  %  is  itself  rate  sensidve,  [54, 
55].  Theoredcal  analyses  of  structural  evoludon  which  can  lead  to  a  rate  sensidve  mechanical 
threshold  were  also  published  recently,  (56, 57], 

An  imponant  observadon,  directly  related  to  ASB  in  XC18  steel,  are  values  of  the  cridcal 
total  strain  (elasuc  plus  plastic)  determined  as  a  funcdon  of  imposed  velocity  on  MDS  specimens. 
The  final  result  is  shown  in  Fig.  3 1  where  the  cridcal  shear  strains.  Fp,  s  5ni/li.  are  shown  as  a 
funcdon  of  log  F.  The  mean  values  are  shown  as  circles.  The  results  for  two  lowest  strain  rates,  f 
•  75  S'*  and  F  •  750  s  ^  were  obtained  with  the  loading  rig  and  the  fast  closed-loop  m^hure,  the 
data  for  three  highest  strain  rates  in  the  range  l.S  ••  10*  s*^  ^  f  ^  4.7  «  10*  s'^  were  obtained  with 
the  new  experimental  technique.  The  impact  velocity  at  the  highest  nominal  strain  rate  was  <- 100 
m/s. 

Since  the  dimensions  of  deforming  parts  of  the  MDS  specimen  are  the  same  as  the  thickness 
of  the  layer  studied  by  numerical  tnethods,  and  the  saoM  steel  was  assumed  for  the  numerical 
study,  a  direct  comparison  is  possible  between  experimental  results  shown  in  Fig.  3 1  and  numeri¬ 
cal  results  of  Fig.  3.  Although  values  localizadon  strains  are  overall  smaller  than  the  critical 
strains  found  by  experiment,  at  lower  strain  rares  the  same  trend  is  observed,  that  is  the  total  shear 
strain  increases  as  a  funcdon  of  log  F.  However,  it  is  clear  that  at  the  highest  strain  rates  an  abrapt 
drop  in  the  cridcal  strain  occurs.  For  example,  the  mean  value  (tf  the  critical  strain  from  six  tests 


perfomed  at  T  »  1.8  *  10^  s*'  (impact  velocity  ~  30  m/s)  is  found  ®  2.60,  whereas  the  mean 
critical  strain  obtained  from  three  tests  at  F  =  4.7  *  10*  s’*  is  Fm  *  0.83.  This  result  clearly  indi¬ 
cates  that  the  mechanism  of  the  .ASB's  formation  and  propagation  differs  substantially  at  different 
initial  and  boundary  conditions.  In  numerical  calculations  boundary  conditions  were  the  same,  and 
were  assumed  as  quasi-steady  state  of  imposed  velocity  with  a  total  absence  of  stress  concentra¬ 
tors.  in  experiment  the  imposed  velocity  rises  from  zero  by  a  projectile  impact  and  stress  concentra¬ 
tors  are  present  as  it  is  shown  in  Fig.  21,  a  and  b.  At  lower  rates  the  experimental  mean  values  arc 
perhaps  higher  in  comparison  to  the  numerical  ones  because  a  small  additional  displacement  has 
occured  due  to  some  plastic  deformation  of  the  specimen  supports. 

A  non-direct  confirmation  of  reduction  of  the  total  shear  strain  at  high  impact  velocity  was 
attempted  by  the  energy  analysis.  The  total  energv-  to  fracture  of  all  MDS  specimens  tested  at  high 
strain  rates  was  detenmined  and  this  result  is  shown  in  Fig.  32.  .Again,  the  energy  drop  is  observed 
at  the  highest  strain  rate  (impact  velocity  -  ICX)  m''s).  The  drop  is  not  too  high  in  comparison  to 
changes  of  the  critical  strain,  mostly  because  of  the  positive  strain  rate  sensitivity  of  flow  stress. 

.A  similar  observation  concerning  energy  consumption  during  fracturing  of  carbon  steels  by 
ASB  formation  was  reponed  in  (58).  Those  results  have  been  re-analysed  by  introduction  of  the 
mean  strain  rate  and  the  tlnai  result  is  shown  in  Fig.  33  for  1010  and  1045  steels  (German  Stan¬ 


dards  CklO  and  Ck43).  The  so  called  "hat”  specimen  geometry  was  employed  and  specimens  were 
deformed  with  SHPB  at  a  narrow  range  of  velocities.  5.5  m-'s  <  V,  s  S.3  m/s.  Again,  it  is  cler 
from  Fig.  33  that  the  fracrur.ng  energy  via  formation  and  propagation  of  .ASB  s  decreases  for  both 
steels  when  the  nominal  strain  rate  is  increased.  If  die  same  loganthmic  scale  were  applied  in  both 
figures  (Fig.  32  and  Fig.  33).  a  drop  of  fractunng  energy  would  be  very  abrupt.  A  similar  beha¬ 


vior.  a  decrease  o:  fracture  smur.  as  .i  function  of  log  f.  was  found  in  {59]  for  93fc  W  alloy.  This 


result  is  shown  in  Fig.  34,  .As  before,  for  .XCIS  steel  the  critical  fracturing  strain  diminishes 
abruptly  around  f  ■  10”*  $'*. 


Tnose  experimental  results  obuuned  very  high  saain  ,'a:e$  are  in  complete  cons^ction  to 
the  numencai  predictions,  for  e.xtunple  Fig.  3.  .A  further  study  \$  .needed  to  clanfy  this  discrepancy. 

Fractographic  study  performei  for  XCl8  steel  by  SE.M  revealed  that  for  all  velocities  the 
fracture  via  .ASB's  occurs  by  creaaon  and  shearing  of  microvoids.  No  so  called  white  bands  wem 
observed  aiTcr  perpendicular  polishing  and  aeaunent  by  .Nital.  Ln  some  areas  a  very  shiny  patches 


ttre  vtstble.  t.hose  are  the  places  where  siidir.g  of  two  neighbouri.ng  pans  of  the  .\IDS  specime.n 


took  place. 

Four  fractograms  for  MDS  specimen  deformed  at  impact  velocity  •  100  m/s  are  shown 
respectively  in  Figs.  35.  36. 37  and  38.  Fig.  35  shows  the  boundary  of  the  shiny  surfKe  and  the 
dimpled  region.  Fig.  36  shows  only  the  shiny  pan  of  the  fracture  surface  with  severe  smeared 
defemnation.  Fig.  37  is  the  dimpled  surface  with  some  traces  of  cleavage.  Finally.  Fig.  38  shows  a 
magnification  of  the  dimpled  region.  Overall,  a  severe  plastic  defomution  occurs  locally  near  dw 
fmcnire  surface  and  some  patches  are  deformed  frirther  due  to  local  sliding.  Temperature  increase 
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in  the  local  sliding  areas  must  be  quite  high,  lliis  is  supponed  by  numerical  results  of  Fig.  3.  The 
local  sliding  and  the  mechanics  of  ASB's  propagation  should  be  studied  further. 

Experiments  with  YAR4340  steel 

In  order  to  study  the  impact  velocity  effects,  or  the  nominal  strain  rate  effects,  on  impact 
shearing,  a  series  of  tests  have  been  performed  with  the  Vacuum  Arc  Remelted  AISI 4340  high 
strength  steel.  The  original  material  was  supplied  by  MTL- Watertown  in  the  form  of  305  x  305  x 
9.52  mm  cross  rolled  plates.  The  chemical  composition  of  the  supplied  VAR4340  steel  was  as 
follows  :  C  0.42  ;  Ni  1.74  ;  Cr  0.89  ;  Mn  0.46 :  Si  0.28  ;  S  0.001 ;  P  0.009  {%  wt.).  About  forty 
.VIDS  specimens  were  machined  with  geometry  shown  in  Fig.  20b.  The  shear  plane  was  assumed 
as  perpendicular  to  the  plane  of  the  original  cross  rolled  piate.  After  precision  machining  the  MDS 
specimens  were  austenirized  1  hour  at  925*C  and  oil  quenched,  re-austeninzed  15  min  at  845®C  and 
oil  quenched  and  finally  tempered  at  200®C  during  2  hours.  The  final  hardness  was  -  52  HRC.  As 
it  has  been  shown  in  [60]  such  thermal  treatment  assures  maximum  of  the  energy  absorbed  up  to 
the  point  of  instability  at  f  ■  10^  S’^ 

Around  30  .VIDS  specimens  were  used  for  quasi-static  and  impact  loading.  In  the  case  of 
loading  at  quasi-stadc  rates  the  rig  shown  in  Fig.  23  was  used  together  with  the  closed-loop  fast 
machine.  In  the  case  of  impact  loading  four  impact  velocities  were  applied :  Vj  ■  20  m/s ;  V2  •  35 
m/s :  V3  •  72  m/s  and  Vj  -  160  m/s.  For  every  nominal  strain  rate  from  three  to  five  specimens 
were  tested.  At  low  strain  rate.  F  «  1.1  *  1(?  s  ’*,  some  specimens  fractured  in  the  mixed  mode, 
that  is  shearing  and  opening.  The  scanning  electron  micrography  (SE.VI)  of  those  surfaces  revealed 
small  facets  of  cleavage  mi.xed  with  some  ductile  tearing.  The  cleavage  facets  showed  a  mean  dia¬ 
meter  of  few  micrometers.  Diameters  of  plastic  dimples  w-ere  much  smaller.  At  all  higher  rates,  and 
also  in  some  cases  at  this  small  strain  rate,  the  hacturing  occured  in  the  shearing  mode. 

The  level  of  shear  strains  to  fracrare  was  quite  limited,  for  example  at  F  ■  l.I  •  10^  s  *  the 
critical  shear  strain  was  around  Fcr  *  0.0425.  The  mean  level  of  the  maximum  shear  soess  for 

.VIDS  specimens  tested  a:  F  «  1. 1  •  10-^  $  ^  was  ■  S63  .MPa.  For  the  same  temper  tested  in 

torsion  (61  ].  KRC  •  54.  at  F  «  6  •  10*^  $4  the  maximum  stress  was  tg,  •  1070  MPa.  It  is  found 
that  relatively  good  agreement  exists  between  data  hom  a  torsion  test  of  a  thin-walkd  ^Kcimen  and 
the  MDS  test  Overall  behavior  of  the  critical  shear  stress  as  a  fiincdoa  of  log  F  is  shown  in  Rgs. 
39a  and  b.  The  mean  values  of  stress  shown  in  Fig.  39b  indkaie  for  a  maximum  of  tc,  at  F  ••  5  • 
IQ)  s'>  and  next  a  drop  at  higher  stnin  rates  occun.  The  same  experimental  results  are  shown  in 
thelinearscakt^FuiFigs.  40a  and  b.  Since  the  maximuffl  impact  velocity  was  Vg*  160  m/s,  and 
at  this  velocity  reladvety  smaH  decrease  of  Tq  found,  it  is  presumed  that  just  beginning  of  a 
more  severe  drop  of  Tct  was  found.  Thus,  the  critical  impact  velocity  may  be  around  200  in/s  with 
the  nominal  velociiy  of  shearing  Fcr  *  ITiose  values  are  hij^ier  in  cotnpanson  to  XC18 


steel.  Fig.  31.  In  ceniin  regions  of  strain  rate  it  is  possible  to  calculate  local  strain  rate  sensitivities 
of  the  critical  stress.  This  is  cemonstrated  in  Fig.  41  and  Fig.  42  where  by  the  broken  lines  the 
linear  regions  are  indicated. 

•At  low  strain  rates  the  standard  rate  sensitivity  may  be  introduced  Pf 


T 


(63) 


The  slope  indicated  by  the  broken  line  gives  Pp  =  54.79  MPa. 

At  high  strain  rate  region  ine  linear  rate  sensitivity  or  pseudo-viscosity  can  be  defined 


the  slope  of  the  broken  line  or  Fig.  42  gives  tIf  =  6.45  Pa.s,  value  companble  to  values  for  XC18 
steel  (Table  21. 

•A  :\-picai  osculogram  for  ‘ne  MDS  specimen  of  4340  steel  tested  at  V,  «=  160  m/s  shown  in 
Fig.  13  indicates  mat  the  total  energy  to  break  the  specimen  is  quite  small  as  compared  to  XCIS 
steel  w  hich  shows  a  Iarge-stn:r.  ?iasrcic>'.  In  the  case  of  .XC18  steel  the  peak  of  the  total  energy  to 
break  the  .MDS  specimen  can  :e  estimated  from  Fig.  32  as  »  230  J.  this  number  gives  me 
energy  absorbed  per  unit  of  sheanng  surface  as  W/.A  «  2.4  J/rrm*.  .An  analysis  of  the  total  energy 
absorbed  for  three  .\IDS  specurer.s  made  of  V.AR43'^  steel  and  loaded  at  impact  velocity  «  160 
m»s.  give  the  mean  value  VV.,A  =  0.50  J/mm-.  The  total  energy  per  unit  area  reponed  in  (63)  to 
break  a  thin-w  ailed  tubular  specimen  of  V.AR4340  steel.  HRC40.  at  F  •  10*  was  W/A  •  0.185 
J.mmv  The  men  e.nergy  acscrcri  up  to  the  point  c:  instabuttv  for  hat"  specimens  of  VAR4340 
steel,  and  for  the  same  temper  as  in  the  present  study,  cxn  be  esumated  from  Fig  7  of  [60j.  At  an 
average  displacement  rate  imposed  on  the  "  hat"  specimen.  V  •  5. 1  ni's.  and  at  the  local  strain  rate 
r  •  10^  .  the  mean  energ>'  per  unit  area  dissipated  up  to  the  rr.a.xiraum  force  is  cstimaicd  as  W/.A 

^  0.084  J/mm-. 

Because  cf  cmferen:  spemoen  geomemKS  and  different  imposed  vesoetties.  as  well  as  thermal 
treaimcRis.  values  of  ^VV.A  :c:  V.AR4340  steel  are  not  equal,  but  they  give  an  order  for  this 
measure.  Also,  utae  to  fra<  ture  were  quite  different  at  different  expehnnenu.  e.xample.  ctse  to 
fracture  in  SHTB  can  be  assumed  as  t^  •  100  ^s.  whereas  time  to  fracture  of  the  MDS  specimen 
can  be  as  short  as  several  iu. 


Since  inioal  and  boundar-  conditions  pUy  an  unporaiu  role  in  ASB's  fonnation  and  fractu- 
n  ig.  comparison  of  rerults  from  different  expenmentai  techniques  must  be  studied  funher. 

The  next  step  in  this  stuc>  was  fractographtc  examination  of  fracrare  surfaces  by  SEM.  .A 
review  of  previous  meullograp.'.tc  and  fractographtc  studies  on  surface  morpherfogies  after  adiaba¬ 
tic  sheanng  was  r.>^n  in  (63\  Ir.  gentral.  there  is  no  agreement  as  to  the  occurence  of  phase 
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transfomiation  in  the  so  called  "white  bands”.  Recent  studies  seem  to  indicate  that  those  bands 
consist  of  a  very  fine-grained,  highly  dislocated  microstructure,  [60].  The  interior  of  so  called 
"transformed"  shear  bands  of  4340  steel  showed  an  extremely  fine  microcrystalline  structure  with 
grain  diameter  ranging  from  8  nm  to  20  nm  [60].  This  observation  supports  the  idea  that  no  trans¬ 
formation  occurs  in  the  white  bands  and  the  interior  of  the  band  is  a  very  heavily  deformed  manen- 
site,  [60].  The  extremely  fine  graines  are  combined  with  the  additional  carbon  in  solution. 

SEM  micrography  should  provide  further  information  on  the  separation  mechanics  of  the 
white  bands.  Figures  43,  44.  45  and  46  show  SEM  pic;  res  of  the  same  region  of  fractured 
surface  of  the  MDS  specimen  deformed  at  impact  velocity  Vo  =  165.8  m/s.  The  three  typical  rrxjr- 
phologies  of  fncture  surfaces  can  be  observed  in  those  figures.  In  Fig.  43,  a  small  magnirication, 
tongues  of  highly  smeared  material  are  identifiable.  The  borders  of  those  tongues  are  marked  by 
the  bnghter  lines.  The  distances  between  the  borders  of  tongues  diminishes  when  the  impact 
velocity  is  increased.  In  Fig.  a  larger  degree  of  magnification,  the  border  of  the  smeared  and 
dimpled  zones  are  clearly  visible.  Fig.  45  shows  the  border  between  smeared  and  dimpled  parts. 
The  end  of  the  tongue  is  covered  with  knobbles,  which  are  characterized  by  spheroidized  surface 
debris,  apparently  from  a  high  temperature,  close  to  the  melting  point  [4:  ].  This  is  probably  a 
separated  white  band  with  a  very  fine  microcrystalline  structure,  this  microscructure  is  partially 
visible  in  the  bottom  of  Fig.  45.  The  area  in  the  upper  side  of  the  tongue,  the  same  magnification  as 
in  Fig.  45.  is  shown  in  Fig.  46.  This  figure  shows  clearly  identifiable  elongated  microvoids  of 
diameter  from  I  M.m  to  3  jim.  The  white  spots  visible  in  most  of  fractograms  are  probably  debris  of 
tefion  dispersed  after  impact. 

.Another  series  of  SE.M  micrograph  examinations  was  perfonned  for  the  MDS  specimen  loa¬ 
ded  with  the  impact  velocitv*  V,  s  71  9  m/s.  Again  three  typical  morphologies  of  fracture  surface 
can  be  identified.  Fig.  47  shows  the  end  of  the  tongue  with  the  white  txvder  line.  The  magnified 
'in  of  the  border  region  is  shown  in  Fig.  4g.  ai  the  left  bottom  the  tongue  is  covered  with 
knobbles,  at  the  right  upper  pan  the  dimples  are  clearly  visible.  The  central  area  of  the  tongue  with 
the  same  magnification  as  the  micropaph  of  Fig.  4$  is  showm  in  Fig.  49.  This  micrograph 
strongly  indicates  for  an  exis-tence  of  a  very  high  temperanire  in  this  region.  The  quasi-dimples  are 
conqsieteiy  spheroidized  and  it  indicates  a  viscous  separation  of  surfaces. 

The  SEM  examination  of  specimen  sur^ces  produced  at  impact  sTlocities  160  m/s  and  72 
trJs  clearly  revealed  patches  of  elonpted  microvoids  (Figs.  46  and  48}.  It  is  presumed  that  auclea- 
liott  and  coakscec^  of  nrerovoids  is  the  {kisary  le^hanism  of  the  surface  sepamtion.  Exissence 
tongues  and  krvjbbled  surfaces  indicatt  for  an  exigence  of  tu^  local  temperaage  (close  to  the 
melting  point)  and  severe  sliding.  Those  rtgions  nuy  be  much  softer,  and  productioa  tongues 
can  be  expUined  by  a  local  dispiaceoient  of  a  soft  maKtiil  by  sUdiag  mfaoes.  When  thidaiess  of  a 
sliding  miocolayer  is  suflic«ndy  large  a  tongue  is  formed  and  die  vicinities  retain  its  dimjrfed  struc¬ 
ture.  Separation  of  tongue  regions  occurs  with  presence  of  very  hi^  temperature,  for  example  Fig. 
49.  Preliminary  analyses  indkaie  for  an  relation  berween  the  mean  tongue  size  and  velocity  of 
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sliding.  Tne  proponion  between  knobbly  and  dimpled  surfaces  may  depend  on  specimen  geome- 
tr>'  and  loading  conditions. 

Present  SEM  fractographic  study  shows  some  similarities  with  results  for  VAR4340  steel 
(HRC40)  reponed  in  [63].  A  further  study  is  needed  to  define  correlanon  between  morphologies  of 
micrographs  and  impact  velocities. 


Discussion  and  conclusions 

This  study  has  raised  some  points  for  funher  studies  in  the  area  of  ASB's  and  fracture 
dynamics. 

There  is  general  agreement  that  the  critical  instability  strain  is  related  to  the  zero  hardening 
condition,  dt/dr  =  0.  of  the  adiabatic  stress-strain  cur\'e.  The  post  instability  behavior  is  much 
more  complex.  In  addition  to  material  behavior,  initial  and  boundary  conditions  could  play  decisive 
role  in  development  and  formation  of  ASB’s  and  also  fracture.  Tne  steady-state  solutions,  like  the 
numerical  one  reported  in  this  report,  or  analytical  ones,  are  very  useful  tools  in  ASB  studies,  they 
can  depict  the  late  stage  morphology  during  ASB  development  as  a  function  of  initial  and  boundary’ 
conditions.  However,  those  solutions  neglect  stress  concentrators  which  are  always  present  in  real 
situations  including  e.xperiment.  It  has  been  shown  in  this  study  that  the  stress  concentrators 
inverse  the  post  iccalization  behavior  as  a  function  of  the  nominal  strain  rate.  To  study  funher 
t.hose  effects  a  variety  of  experimental  techniques  must  be  applied  with  well  specified  initial  and 
boundary  conditions.  The  next  inevitable  step  is  a  superposition  of  fracture  dynamics  on  ASB 
morphology.  Comparison  of  all  available  approaches  for  the  same  or  similar  materials  may  provide 
a  data  base  on  development  and  formation  of  ASB’s, 

It  may  be  argued  that  in  order  to  compare  the  susceptibility  to  adiabatic  shear  in  different 
conditions,  that  is  in  presence  of  different  stress  concentrators,  it  is  necessary  to  define  the  defect 
that  initiates  the  localised  .ASB.  It  was  suggested  in  [64]  that  the  defect,  or  stress  concentrator, 
must  be  scaled  in  an  appropriate  way  to  the  homogeneous  response  in  the  onset  of  ASB. 

The  new  experimental  technique  developed  during  this  project  permits  for  shear  testing  at 
large  strains  within  a  wide  range  of  the  nominal  strain  rates,  from  quasi-static  up  to  5  x  10^  s*^ 
Tne  local  strain  rates  during  localization  must  be  much  higher.  The  configuration  of  impact  experi¬ 
ment  with  MDS  specimen  and  direct  optical  measurement  of  displacement  permits  for  sufficient 
flexibility  to  change  the  notch  geometry.  The  limiting  case  is  to  apply  MDS  specimen  with  symme¬ 
tric  fatigue  cracks  to  study  Mode  n  fracture  under  impact  conditions.  A  general  study  of  Mode  II, 
quasi-static  loading  and  dynamics,  has  been  recently  completed  as  the  Appendix  N®  1  to  this 
Repon,  [4],  Dynamic  plasticity  has  been  included  in  this  Mode  II  study.  Also  differen*  notch 
geometries  could  be  applied  to  vary  stress  concentrators. 
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The  most  imponant  finding  of  this  study  is  confirmation  of  the  existence  of  the  critical  impact 
velocity  in  shear  deformation.  This  velocity  is  around  100  m/s  for  the  mild  steel  XC18  and  above 
160  m/s  for  VAR4340  steel  (HRC  52).  The  experimental  technique  with  MDS  specimens  permits 
for  still  higher  impact  velocities  to  be  tested. 

Another  imponant  area  in  studying  ASB's  is  post  factum  metallographic  examination  of  frac¬ 
ture  surfaces.  In  the  present  study  the  SEM  has  been  applied  for  both  steels  studied  (XC18  and 
VAR4340  -  HRC52).  The  mild  steel  XC18  showed  large  stri  ins  before  the  onset  of  instability.  But 
those  large  strains  were  substantially  reduced  at  impact  velocity  ~  100  m/s.  In  the  case  of  this  steel 
a  substantial  adiabatic  increase  of  temperature  occurs  not  only  during  localization  but  before.  On 
the  contrary,  in  the  case  of  VAR4340  the  temperature  increase  prior  to  the  onset  of  instability  is  not 
high,  but  during  an  abrupt  localization  the  melting  point  can  be  reached  in  rubing  islands,  Fig.  49. 
The  morphology  of  fracture  and  shearing  velocity  should  be  studied  further,  using  perhaps  an 
image  analysis. 

Although  considerr.bly  more  evidence  must  be  gathered  for  other  materials  and  other  defect 
geometries,  this  study  has  revealed  that  the  very  imponant  parameter  in  studying  ASB's  is  the 
impact  velocity. 
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Figures 

Fig.  1  a.  Gradient  of  temperature  OT/8x)  as  a  function  of  log  for  L  =  2.0  mm  ;  1-miId  steel, 
2-aiuminum,  3-copper ; 

b.  Values  of  c  as  a  function  of  log  Fn  ;  L  =  2.0  m,  maximum  of  \  defines  transidon 

(Fn)ir- 

Fig.  2  Criricai  values  of  shear  strain  rates  (rn)cr  as  a  funedon  of  tube  length  L  ;  +  -  copper, 
*  -  aluminum,  o  -  steel. 

Fig.  3  Nominal  values  of  instability  strain  Fc  and  localizadon  strain  Floc  as  a  funedon  of  the 
nominal  strain  rate  F;,  =  V/h,  numerical  simulation  for  XC18  steel. 

Fig.  4  Nominal  values  of  instability  strain  Fc  in  expanded  scale  as  a  funedon  of  the  nominal 
strain  rate  F-.  numencai  simulation  for  XC18  steel. 

Fig.  5  The  rado  of  F^F -s  a  function  of  log  F^.  numerical  simulation  for  XC 1 8  steel. 

Fig.  6  Localization  of  shear  strain  at  different  levels  of  the  nominal  strain  F-  at  F^  =  I(P  s*L 

Fig.  7  Localizadon  of  temperature  at  different  levels  of  the  nominal  strain  Fn  at  Fa  =  10^  s‘^ 

Fig.  8  Evolution  of  srain  gradients  gr  (y)  at  the  same  levels  of  the  nominal  strain  as  indicated  in 
Fig.  6.  nomtnai  strain  rate  f*  =  10^  s*\ 

Fig.  9  Evolution  of  strain  gradients  in  the  halt-thickness  of  the  layer  at  F-  =  10^  s'‘  and  diffe¬ 
rent  icvcis  of  nominal  scrain  F- ;  dashed  line  shows  me  asymptotic  procedure  in  finding 
the  ASB  thic  .ness 

Fig.  10  Evolution  of  the  .ASB  band  thickness  as  a  funedon  of  me  nominal  srain  rate  ;  o  -  analysis 
based  on  srain  gradiena  :  •  -  analysis  based  on  temperature  gradients. 

.New  configuration  of  e.vpehmental  setup  developed  :n  LPMM-Metz  ;  P  -  projectile  : 
S  -  double  "shea:  spcc;mc.n  .  L  -  source  of  light ;  1.2.3  -  fiber  optics ;  F  -  photodiodes : 
TCI.  TC2  •  time  counters .  E  •  optical  extensometer ;  Tr  •  strain  resistance  gage  :  S  -  DC 
supply  unit :  .Ai  •  ampiificr ;  DO  •  digital  osciiioscope  :  PC  •  personal  computer ;  R  •  XT 
recorcer  or  gnph:c  pr.nter. 

Oscillogram  from  me  ast  of  XCIS  sreel  at  Lmpact  velocity  a  35.0S  nv’s ;  upper  trace  is 
tnc  signal  from.  SR  station  TI  :  lower  trace  is  the  dispiaccmcr.;  signal  from  the  opdcil 
extensometer. 

Fig.  1 3  Oscillogram  from  the  test  of  VaR  4340  steel  at  impact  velocity  •  160  nvs  ;  the  upper 
trace  is  die  st^al  from  SR  stauon  Tl.  the  low-er  trace  is  the  displacement  signal  from  the 
opt^  extensometer. 

Fig.  14  Deaiis  of  the  Hopkinson  tube ; 

a.  the  face  at  which  pressure  is  applied  •  shadowed  surface  supports  MDS  specimen 
ia  X  bi : 

b.  dimensions  ami  disunces  used  in  FE  calculaoons  x.  »  90.5  mm ;  x:  s  170.5  mm. 


r:2. 
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Fig.  15  Initial  conditions  of  pressure  on  the  surface  2  (a  x  b)  as  sliown  in  Fig.  14  ;  trapezoid 
pulses  with  different  times  tr,  to  and  t,  as  shown  at  the  bottom. 

Fig.  16  Results  of  FE  calculations  e  (t)  for  xi  =  90.5  mm,  case  N®  1  of  Fig.  14,  solid  and  bre 
ken  lines  are  respectively  for  external  and  internal  surface  of  the  P»be ;  a  -  plane  parallel  tc 
the  MDS  specimen ;  b  -  perpendicular  plane. 

Fig.  17  Results  of  FE  calculations  e  (t)  for  X2  =  17C.5  mm,  case  N®  1  of  Fig.  14,  solid  and  bro¬ 
ken  lines  are  respectively  for  external  and  internal  surtace  of  the  mbe ;  a  -  plane  parallel  to 
the  MDS  specimen ;  b  -  perpendicular  case. 

Fig.  18  Results  of  FE  calculations  £  (t)  for  x\  =  90.5  mm,  :.ase  N®  4  of  Fig.  14. 

Fig.  19  Results  of  FE  calculations  £  (t)  for  X2  =  170.5  mm,  case  N®  4  of  Fig.  14. 

Fig.  20  Two  specimens  with  modified  geometries  for  double  shear  tests ;  a  -  geometry  initially 
assumed  after  modifications  with  respect  tc  [48] ;  b  -  perfected  geometry  with  enlarged 
supports. 

Fig.  21  Distribution  of  shear  stresses  in  MDS  sp'cimen  after  nominal  shear  strain  ;  (a)  F  =  2.5  x 
10-3 ;  (b)  r  =  0.50. 

Fig.  22  Distribution  of  the  shear  strain  component  in  MDS  specimen  af'ter  nominal  shear  strain  : 
(a)  r  =  2.3  X  10-® ;  (b)  r  =  0.50. 

Fig.  23  Rig  for  deformation  of  .MDS  specimens  at  lifferent  strain  rates,  from  10"^  s-‘  to  lO’ s  *, 
using  a  fast  closed-loop  testin"  machine ;  MDS  specimen  is  shadowed. 

Fig.  24  Values  of  the  nominal  shear  strain  rule  F  as  a  function  of  gage  length  U  for  two  a.xial 
velocities  20  m/s  and  i  <0  m  s ;  CF  •  gage  lengths  used  in  [48] ;  JRK  gage  length  1, 
=  2.0  mm  used  in  MDS  specimen. 

Fig.  25  Variation  of  upper  yield  stress  t»  in  shear  with  logarithm  of  strain  rate,  log  F,  MDS  spe¬ 
cimen. 

Fig.  26  Upper  vicld  stress  ty  versus  shear  itrain  rate  F  in  linear  scale  of  F. 

Fig.  27  Variation  of  lo>ker  yield  stress  tyj  in  shear  with  logarithm  of  shear  strain  rate,  log  f ;  •  - 
MDS  specimen :  •  experimental  points  horn  ref.  (48). 

Fig.  28  Comparison  of  results  for  lower  yield  stress  tyi  for  XC18  steel  of  Modified  Double  shear 
(>1,  and  inclined  impact  (o;  for  pure  iron.  ( ISj. 

Fig.  29  Upper  yield  stress  Xy  versus  shear  strain  rate  in  linear  scale,  XC18  steel.  MDS  sped- 
men ;  broken  line  shows  the  best  liirear  fit 

Fig.  30  Lower  yield  stress  tyi  versus  shear  strain  rate  T  in  linear  scale :  XC18  steel.  MDS  speci¬ 
men  ;  broken  line  shows  the  best  litKar  fit 

Fig.  31  Shear  strains  of  localization  vent‘d  log  F :  XC18  steel.  MDS  specimen  ;  circles  denote 
mean  values. 
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Fig.  32  Tocal  fracturing  energy  for  MDS  specimens  as  a  function  of  log  f ,  range  of  high  strain 
rates,  steel  XC18. 

Fig.  33  Fracturing  energy  for  "hat"  speciments  as  a  function  of  log  F,  results  from  [58]  re-analy¬ 
sed  for  two  carbon  steels. 

Fig.  34  Fracture  strain  as  a  function  of  strain  rate  for  93%  W  alloy,  after  [59]. 

Fig.  35  SEM  micrograph  of  transition  region  from  smeared  shiny  surface  to  microvoid  separa¬ 
tion,  steel  XC18,  Vo  *  100  m/s. 

Fig.  36  SEM  micrograph  of  smeared  shiny  surface,  steel  XC18,  Vq  *  100  m/s. 

Fig.  37  SEM  micrograph  of  region  of  heavilv  elongated  dimples  with  traces  of  cleavage,  steel 
XC18,  lOOm/s.  " 

Fig.  38  SEM  micrograph  showing  magnified  pan  of  Fig.  37,  steel  XC18,  V,, »  100  m/s. 

Fig.  39  Cridcal  shear  stresses  t;;  for  V..\R  4340  steel  versus  log  F ;  a  -  data  from  inoividual  MDS 
tests ;  b  -  mean  vaiues. 

Fig.  40  Critical  shear  stresses  tor  V.\R  4340  steel  versus  linear  scale  of  f ;  a  -  data  from  indi- 
\nduai  tests ;  b  •  mean  vaiues. 

Fig.  41  Results  for  V.Ml  4340  steel  indicating  appro.ximadon  by  rate  sensitivity  (Jp,  eq.  (63). 

Fig.  42  Rssuits  tor  VAR  4340  steel  showing  approximation  by  pseudo-nscosity  rjp,  eq.  (641. 

Fig.  43  SEM  nticrograph  of  fracrore  sunace.  .MDS  specimen  Vq  =  165.8  nvs. 

Fig.  44  SE.M  micr'^graph  of  fracture  surface  with  higher  magnification,  V'o  =  165.8  tn/s. 

Fig.  45  ShM  micrograph  of  fracture  surface,  end  of  tongue.  V,  s  165.8  m/s. 

Fig.  46  SE.M  micrograph  of  fractu.’e  surface,  dimpled  region.  =  165.8  m  s. 

Fig.  4“  SE.M  micrograph  of  fracrare  su.nace.  region  of  tongue,  impact  velocity  =  71.9  m/s. 

Fig.  4S  SE.VI  micrograph  of  fracture  sunace.  end  of  tongue,  =  71.9  nvs. 

Fig.  49  SEM  micrograph  of  fracrure  sunace.  central  pun  of  die  tongue  in  Fig.  47.  =  T 1 .9  m/s. 
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